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Abstract 

This  study  investigated  seasonally  demanded  consumable 
items  at  the  base-level.  This  study  examined  how  the  Stan¬ 
dard  Base  Supply  System  currently  addresses  seasonally 
demanded  consumable  items  and  some  alternative  methods  of 
addressing  consumable  seasonal  demand  items  in  the  SBSS . 

The  alternative  methods  analyzed  in  depth  were  simple  and 
Winters'  seasonal  exponential  smoothing,  and  Box-Jenkins 
forecasting  models. 

This  study  found  that  items  under  study  display  some 
seasonal  demand  tendencies.  The  sample  consisted  of  12  out 
of  77  items  identified  as  seasonal  by  Stock  Control  person¬ 
nel  at  Langley  AFB ,  VA .  A  graphical  analysis  showed  strong 
er  seasonal  demand  tendencies  than  did  the  autocorrelation 
function  in  which  the  correlations  between  demands  one  year 
apart  are  determined. 

As  was  expected,  the  two  seasonal  models.  Winters' 
exponential  smoothing  and  Box-Jenkins  better  predicted 
demands  for  items  under  study  than  the  SBSS  model.  Of  28 
forecasts.  Winters'  exponential  smoothing  was  best  13  times 
while  Box-Jenkins  models  were  best  9  times. 

The  autocorrelation  function  could  be  used  to  test 
demand  data  for  seasonality  and  flag  items  with  seasonal 
demand  patterns  for  special  seasonal  treatment,  but  this  is 
not  currently  practical.  Any  useful  effort  to  test  all 


xi 


items  loaded  on  the  SBSS  at  a  base  would  require  demand  data 
for  each  item  for  four  or,  preferably,  more  years.  •  Any 
manual  attempts  to  test  items  for  seasonality  would  be 
impractical  given  the  number  of  items  in  the  average  base 
supply  account . 

Among  the  recommendations  given  as  a  result  of  this 
study  is  the  suggestion  that  additional  work  be  done  to 
facilitate  basing  demand  forecasts  on  seasonal  models  where 
appropriate . 
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THE  IDENTIFICATION  AND 


FORECASTING  OF  SEASONAL  DEMAND 
CONSUMABLE  ITEMS  IN  BASE  SUPPLY 


I .  Introduction 


General  Issue 

In  1985,  Headquarters  Tactical  Air  Command  Supply  pei — 
sonnel  identified  to  the  Air  Force  Logistics  Management 
Center  the  need  "to  develop  an  automated  method  to  identify 
seasonal  items  "(1).  In  the  Standard  Base  Supply  System 
(SBSS) .  there  is  no  method  to  identify  consumable  items  with 
seasonal  demand  patterns  (1).  This  inability  to  determine 
periods  of  high  and  low  demands  leads  to  excessive  inven¬ 
tories  of  items  during  periods  of  low  demand  and  stockouts 
during  periods  of  high  demand  (1).  While  a  method  exists  in 
the  supply  system  to  load  special  stockage  levels  for  sea¬ 
sonal  items,  no  automated  system  is  available  to  identify 
those  items  or  the  appropriate  special  levels  (1). 

Specific  Problem 

At  the  base-level,  managers  are  forced  to  use  guesswork 
and  personalized  attention  to  forecast  demand  and  ensure 
availability  of  seasonally  demanded  items.  In  this  context, 
seasonal  items  are  considered  items  with  demand  patterns  of 
three  or  more  consecutive  months’  demand  exceeding  one 
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standard  deviation  from  the  mean  demand.  Once  an  item  is 
identified  as  seasonal,  managers  must  also  determine  the  ap¬ 
propriate  stockage  level  for  the  item. 

Research  Questions 

The  following  areas  were  addressed  to  determine  if  an 
automated  system  of  forecasting  range  and  depth  of  consum¬ 
able  seasonal  demand  items  (what  items  to  handle  as  seasonal 
and  how  many  of  those  items  to  stock)  can  be  developed  for 
the  SBSS : 

1.  How  does  the  SBSS  currently  address  consumable  sea¬ 
sonal  demand  items? 

2.  What  are  alternative  methods  of  addressing  consum¬ 
able  seasonal  demand  items  in  the  SBSS? 

3.  Do  items  currently  identified  as  seasonal  by  person¬ 
nel  at  the  base-level  actually  display  seasonal  tendencies? 

4.  Can  demand  for  the  items  currently  identified  as 
seasonal  by  personnel  at  the  base-level  be  better  forecast 
using  a  model  that  incorporates  seasonal  demand  pattern 
information,  or  the  current  SBSS  model? 

5.  What  method  could  be  used  by  personnel  at  the  base- 
level  to  identify  items  as  seasonal  and  forecast  demands  for 
those  items  accordingly? 

Scope 

This  research  effort  concentrates  on  consumable  items 
(also  known  as  consumption  items) 
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which  are  either  consumed  in  use  or  which  lose 
their  original  identity  during  periods  of  use  by 
incorporation  into  or  attachments  upon  another 
assembly.  Consumption  items  are  issued  on  an  as 
required  basis  and  consist  of  such  supplies  as 
maintenance  parts,  raw  materiel,  office  or  house¬ 
keeping  supplies  consumed  in  use,  and  other  sim¬ 
ilar  items.  (9:158-159) 

Consumable  items  generally  are  more  economical  to  replace 
than  to  repair.  (Parkas,  non-consumable  items,  are  in¬ 
cluded.  in  this  study  because  they  were  included  in  the 
sample  of  seasonal  items  by  the  test  base.) 

Since  the  SBSS  maintains  limited  demand  data,  the  re¬ 
search  is  limited  to  one  base  which  contributes  to  the  Air 
Force  Logistics  Management  Center's  (AFLMC)  database.  Ini¬ 
tially  three  of  the  twelve  contributing  bases  were  selected 
according  to  several  criteria:  1)  their  location  in  an  area 
with  four  distinct  seasons,  2)  being  in  different  major 
commands,  and  3)  participating  in  the  AFLMC  database.  flic 
bases  selected  as  potential  candidates  were  Minot  Air  Force 
Base,  ND,  of  the  Strategic  Air  Command,  Dover  Air  Force 
Base,  DE,  of  the  Military  Airlift  Command,  and  Langley  Air 
Force  Base,  VA ,  of  the  Tactical  Air  Command.  Because  of  the 
effort  involved  in  extracting  data  from  the  AFLMC  database 
it  was  decided  to  limit  the  study  to  only  one  base,  Langley 
AFB .  VA. 


Summary  and  Overview 

This  chapter  presented  the  problem  of  a  lack  of  a  sys¬ 
tematic  method  for  identifying  seasonal  items  in  the  base 


supply  system.  The  chapter  also  described  the  specific 
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problem,  reviewed  the  research  questions,  and  delineated  the 
scope  of  the  research.  The  next  chapter  will  cover  the 
requirements  computation  and  demand  forecasting  processes  of 
the  SBSS .  In  addition,  a  brief  overview  of  several 
forecasting  methods  will  be  presented.  The  second  chapter 
concludes  with  a  survey  of  previous  studies  in  this  area. 
Chapter  Three  continues  by  detailing  the  research  methodol¬ 
ogy.  Chapter  Four  presents  the  results  and  analysis  of  the 
data  collected.  Finally,  Chapter  Five  gives  the  conclusions 
and  recommendations  derived  from  the  research. 
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II.  Literature  Review 


Introduct ion 


The  United  States  Air  Force  is  one  of  the  largest  buy¬ 
ers  of  goods  and  materials  in  the  world.  In  consumable 
items  alone,  over  519.000  different  items  are  stocked  in  the 
AFLC  inventory  for  support  of  weapons  systems  (10).  With 
such  large  purchasing  needs  lies  the  inherent  responsibility 
to  manage  assets  in  an  "effective  and  efficient"  manner 
(8:2-45) . 

Once  an  item  is  purchased,  it  is  either  used  or  held  in 
inventory  until  needed  or  deemed  excess  to  requirements.  As 
Ammer  notes,  inventories  act  as  a  hedge  against  uncertain¬ 
ties  in  supply  and  demand  (2:257).  He  also  notes  that  with 

too  little  inventory  .  .  .  manufacturing  [or  pur¬ 

chasing]  efficiency  and  customer  relations  are 
bound  to  be  hurt.  Stockouts  of  essential  mate¬ 
rials  mean  some  interruption  of  production  [or 
readiness],  which  raises  costs  .  [In  con¬ 

trast,  too  much  inventory  will]  tie  up  a  company's 
capital:  they  generate  storage  costs;  and  deterio¬ 
rate  or  become  obsolete  in  storage.  (2:255-256) 


Standard  Base  Supply  System  Requirements  Computations 

At  the  base,  or  retail,  level  the  SBSS  is  a  computer¬ 
ized  system  designed  to  monitor  inventories  of  spare  parts, 
equipment  and  goods  (8:1-5).  According  to  Patterson.  Sep- 
panen  described  the  SBSS  "as  a  multi-item,  single  echelon, 
continuous  review  inventory  system  with  stochastic,  multiple 
unit  demands,  backordering  and  an  annual  budget  constraint" 
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(23:1).  This  means  the  system  manages  a  multitude  of  items 
at  a  single  (base)  level,  and  has  continuous  knowledge  of 
stock  levels.  The  system  also  assumes  demand  is  stoch¬ 
astic.  It  permits  backorders,  and  allows  demands  for  multi¬ 
ple  units,  all  subject  to  a  budget  constraint.  Patterson 
also  observed  that  the  daily  demand  rate  (DDR)  "is  the  fore¬ 
cast  measurement  for  the  SBSS .  Unless  accurate  estimates 
are  available  for  the  DDR  considerable  errors  may  result" 
(23:3) .  Thus,  a  critical  function  of  the  base  supply  system 
is  to  forecast  demands  to  determine  the  quantity  of  each 
item  to  hold  at  the  base-level.  The  DDR,  or  estimated 
demand,  is  the  backbone  of  the  entire  SBSS  stockage  policy. 

For  years,  logisticians,  both  civilian  and  military, 
have  been  trying  to  predict  demand  based  on  historical  data. 
As  early  as  1962.  Solomon  noted  that  demand  for  Navy  air¬ 
craft  parts  was  exceedingly  "low  and  sporadic,"  making 
forecasting  difficult  (28:55).  In  light  of  this,  the  SBSS' 
method  of  determining  forecasts  is  discussed,  along  with  how 
those  forecasts  are  used.  A  brief  survey  of  other  fore¬ 
casting  methods  will  follow.  Finally,  a  close  examination 
of  two  powerful  forecasting  methods,  exponential  smoothing 
and  the  Box-Jenkins  methodology  will  ensue. 

Demand  Forecasting  in  the  SBSS.  As  was  mentioned  pre¬ 
viously.  demand  computations  are  a  critical  component  in  the 
calculation  of  stockage  requirements.  The  DDR  is  used  to 
determine  both  the  quantity  to  stock  and  which  items  to 
stock  (23:1).  First  the  calculations  involved  in 
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determining  the  quantity  to  stock  will  be  presented.  The 
USAF  Supply  Manual  states  that  consumable  item  demand  levels 
"are  based  on  an  economic  order  policy  that  balances  the 
cost  to  order  with  cost  to  hold"  (8:19-15).  For  a  given 
item,  the  demand  level  is  calculated  from  the  sum  of  the 
"economic  order  quantity  (EOQ).  order  and  shipping  time 
(0&3T)  quantity,  safety  level  quantity  (SLQ) ,  and  a  0.999 
rounding  factor"  (8:19-15).  The  demand  levels  are  updated 
quarterly,  for  stocked  items,  in  order  to  compensate  for  any 
changes  in  observed  demand  (8:19-11). 

Daily  Demand  Rate.  Common  to  the  calculation  of  the 
EOQ,  O&ST  quantity,  and  SLQ  is  the  daily  demand  rate  (DDR). 
The  daily  demand  rate  formula  is: 

DDR  -  Cumulative  Recurring  Demands _  [ 1 ] 

Current  Date  -  Date  of  First  Demand  (8:19-51) 

The  actual  forecasting  technique  is  more  complex  than 
first  appears  from  this  formula,  since  both  the  Date  of 
First  Demand  (DOFD)  and  Cumulative  Recurring  Demands  (CRD) 
are  updated  periodically.  The  CRD  is  simply  a  counter  that 
records  the  number  of  total  demands  for  the  item.  It  is 
increased  by  the  number  of  units  ordered  for  each  demand 
(23:2).  In  addition,  the  denominator  in  [1]  is  replaced 
with  180  (six  months  in  days)  when  less  than  180  days  of 
demand  history  (current  date  -  DOFD)  is  available.  This 
assures  that  new  items  are  not  overstocked  (23:2).  Two 
additional  adjustments  are  made  every  six  months.  The  CRD 
is  updated  to  the  current  DDR  times  the  less.er  of  the 
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difference  between  the  current  date  and  DOFD,  or  365  (23:2). 
Also,  the  DOFD  is  adjusted  to  be  the  maximum  of  the  DOFD  or 
the  difference  between  the  current  date  and  365  (23:2).  The 
net  result  of  these  adjustments  is  that  the  item's  DDR  is 
based  on  up  to  540  days  of  demand  history,  (23:3)  and  that 
the  forecasting  model  becomes  "a  modified  exponential 
smoothing  (model!  with  a  changing  smoothing  parameter" 
(23:3).  The  exponential  smoothing  model  is  presented  in 
greater  detail  later  in  this  paper. 

Economic  Order  Quantity.  The  economic  order  quantity 
is  a  model  that  replaces  stock  in  batch  sized  orders  to 
minimize  the  sum  of  carrying  cost  (the  cost  to  hold  an  item 
in  inventory)  and  the  ordering  cost.  The  simple  EOQ  model 
requires  several  basic  assumptions  including  (20:638): 

1)  The  demand  rate  is  fixed. 

2)  Backorders  or  shortages  are  not  allowed. 

3)  Delivery  occurs  just  as  inventories  are  depleted. 

4)  The  quantity  ordered  does  not  affect  carrying  and 
holding  costs. 

5)  Item  cost  is  fixed. 

6)  Orders  of  one  item  do  not  affect  another  item's 
ordering  or  holding  costs. 

The  simple  EOQ  formula  then,  is: 

Q-(2DC  /Ch)H  [21 

(20:643) 

where  Q  is  the  optimal  order  quantity,  D  is  the  annual 
demand,  CQ  is  the  cost  to  place  an  order,  and  is  the  cost 
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to  hold  one  unit  of  the  item  in  stock  for  one  year.  Holding 
cost.  Ch  is  usually  expressed  as  a  percentage  of  the  cost  of 
the  item  (20:643).  So.  the  EOQ  is  based  on  four  factors: 
the  cost  to  order  the  item,  the  cost  to  hold  the  item  in 
inventory,  the  annual  demand  for  the  item,  and  the  purchase 
price  of  the  item  (8:19-15). 

Air  Force  SBSS  stockage  policy  sets  the  cost  to  order 
at  two  values  depending  on  the  source  of  the  item.  For 
items  centrally  procured  from  the  depot  supply  system  the 
cost  to  order  is  set  at  $5.20  (8:19-15).  For  local  pur¬ 
chase  items  it  is  $19.94  because  of  the  added  costs  of 
contracting  to  buy  the  item  on  the  local  market  (8:19-15). 
The  cost  to  hold  one  unit  in  stock  for  one  year  is  set  at 
15%  of  the  item's  purchase  cost  (8:19-15). 

This  leads  to  the  actual  SBSS  EOQ  formula.  Since  the 
cost  to  hold  and  cost  to  order  are  assumed  constant  for  both 
purchase  cost  categories,  two  formulas  are  used,  one  for 
centrally  procured  items  and  one  for  local  purchase  items: 


EOQ(Nonlocal  Purchase)  =  , 

[8.3  *(DDR  *  365*  Unit  Price)"5  ]/Unit  Price  [3] 

(8:19-52) 

EOQ (Local  Purchase)  « 

[16.3  *(DDR  *  365  *  Unit  Priced  ]/Unit  Price  [4] 

(8: 19-52) 

The  EOQ  is  not  used  in  the  strictest  sense,  however,  since 
Air  Force  stockage  policy  dictates  that  the  EOQ  value  must 
fall  between  30  times  the  DDR  and  365  times  the  DDR.  (8.19- 
52).  So,  at  times,  the  exact  EOQ  quantity  may  not  be 
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ordered.  In  addition  to  the  EOQ  formula,  however,  several 
other  quantities  must  be  computed  before  determining  the 
total  quantity  to  stock  of  an  item. 

Order  and  Ship  Time  Quantity.  The  first  of  these  other 
quantities  is  the  order  and  ship  time  quantity.  This  is 
simply  the  product  of  the  DDR  and  the  average  shipping  time 
in  days  between  the  source  of  supply  and  the  base  (8:19- 
53) .  This  value  gives  the  quantity  to  stock  to  compensate 
for  demands  during  shipment  time  from  the  depot. 

Safety  Level  Quantity.  The  second  quantity  is  the 
safety  level,  or  the  amount  of  stock  held  to  compensate  for 
variations  in  both  lead  time  demand  and  order  and  shipping 
times.  First  the  statistical  variance  of  demand  (VOD)  and 
the  variance  of  order  and  shipping  time  (V00)  must  be  calcu¬ 
lated.  To  understand  the  following  formulas,  and  their 
terms,  several  ideas  must  be  introduced.  The  VOD  is 
straightforward:  it  is  simply  a  measure  of  how  much  lead 

time  demand  may  vary  over  time  for  the  item.  The  VO 0. 
however,  is  calculated  for  all  items  coming  from  a  given 
source  of  supply.  For  each  receipt  of  an  item  from  the 
source  of  supply,  the  time  it  took  to  receive  the  item  is 
recorded.  For  example,  if  it  took  an  item  12  days  to  ar¬ 
rive,  then  the  counter  for  arrivals  in  1  to  15  days  (a 
delivery  time  segment)  for  that  source  of  supply  would  be 
increased  by  one.  Then  a  reca 1 culat ion  according  to  the 
following  formulas  would  occur: 
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VOD  -  2 (Demand^)-  (2Demand) ^ 


n1  (8:19-53) 


VOO  =  2FI*MI2  -  ( 2FI *MI ) 2 

n0 

- b -  [6] 

n2  (8:19-53) 

where  Is  the  maximum  of  either  the  number  of  days  from 
the  first  demand  or  180.  n->  is  the  number  of  receipts  of  the 
item  for  the  period  under  consideration,  FI  is  the  number  of 
receipts  for  a  particular  delivery  time  segment  coming  from 
a  given  source  of  supply,  and  MI  is  the  midpoint  of  the 
delivery  time  segment  days  allowed  for  shipping  for  a  par¬ 
ticular  priority  coming  from  a  given  source  (8  for  the  time 
segment  1-15  days)  (8:19-53).  The  VOD  is  calculated  during 
the  requirements  computation  each  time  an  item  is  demanded 
or  quarterly  during  file  status.  The  VOO  is  calculated 
quarterly  for  items  with  100  or  more  receipts.  These  values 
are  then  used  in  the  formula 

SLQ  =  C  *  (O&ST  *  VOD  +  DDR2*  VOO) 1/2  [7] 

(8: 19-53) 

where  C  is  a  multiplier  which  increases  the  safety  level  of 
stocks  according  to  the  perceived  criticality  of  the  unit's 
mission  and  SLQ  is  the  safety  level  quantity,  the  amount  of 
stock  to  hold  to  prevent  stockouts  due  to  variations  in 
shipping  time  and  lead  time  demand.  (O&ST  as  used  in  this 
equation  is  the  mean  of  historical  O&ST ' s  recorded  for  items 
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with  a  particular  priority  coming  from  a  given  source  of 
supply . ) 

Economic  Order  Quantity  Demand  Level.  Once  the  neces¬ 
sary  parameters  are  calculated,  the  economic  order  quantity 
demand  level  may  be  calculated  as  follows: 

EOQDL=  TRUNC  [EOQ  +  O&STQ  +  SLQ  +  .999]  [Q] 

(8: 19-52) 

where  TRUNC  means  only  the  integer  portion  of  the  number  is 
retained  and  the  .999  serves  as  a  round  up  factor.  The 
EOQDL  may  be  considered  the  amount  of  stock  to  hold  which 
minimises  cost  of  holding  and  ordering  while  limiting  stock¬ 
outs  due  to  fluctuations  in  shipping  time  and  lead-time 
demand . 

Thus,  as  equations  3,  4.  5,  and  7  show,  the  daily  de¬ 
mand  rate  (estimated  demand)  is  the  backbone  of  the  entire 
SBSS  stockage  policy.  Blazer  noted  that  "The  forecast  for 
demand  is  the  biggest  factor  in  determining  the  depth  of 
stock"  (4:1).  It  appears  in  the  basic  EOQ  formula,  the 
order  and  ship  time  calculation,  and  the  safety  level 
computation.  Substantial  errors  in  this  value  would 
certainly  lead  to  potentially  costly  errors  in  stockage 
levels,  either  in  terms  of  monetary  expenditures  or  mission 
support . 
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Forecasting  Methods 

Because  of  the  importance  of  the  demand  forecast,  the 

methods  by  which  this  value  can  be  determined  must  be 

examined.  A  number  of  methods  of  forecasting  demand  are 

available.  Fildes  divides  them  into  three  major  categories: 

The  judgmental — where  individual  opinions  are  pro¬ 
cessed,  perhaps  in  a  complicated  fashion. 

The  extrapolative — where  forecasts  are  made  for  a 
particular  variable  using  only  that  variable's  past 
history.  The  patterns  identified  in  the  past  are 
assumed  to  hold  over  to  the  future [.] 

The  casual  (or  structural) — where  an  attempt  is 
made  to  identify  relationships  between  variables 
that  have  held  in  the  past,  for  example,  volume  of 
brand  sales  and  that  product's  relative  price.  The 
relationships  are  then  assumed  to  hold  into  the 
future .  (11:92) 

Judgmental  Forecasts.  Fildes  divides  judgmental  fore¬ 
casts  into  three  categories:  individual,  committee  and 
Delphi  forecasts  (11:93).  Individual  forecasts  are  made 
when  a  single  person  forms  an  opinion  on  future  demand  based 
on  experience  and  familiarity  with  the  current  situation 
(11:93).  Committee  or  survey  forecasts  may  be  thought  of  as 
occurring  when  a  group  of  people,  either  management  or 
customers,  form  an  opinion  on  future  demand  based  on  their 
perceptions  (11:93).  The  Delphi  method  uses  special  tech¬ 
niques  to  combine  the  ideas  of  several  people,  but  gives  the 
individuals  feedback  without  allowing  one  person's  ideas  to 
dominate  the  decision  process  (11:93). 

Extrapolative  Forecasts.  According  to  Fildes,  extrapo¬ 
lative  forecasts  come  in  five  basic  forms  (11:94).  The 
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first  is  trend  curve  analysis,  the  method  of  fitting  past 
demand  data  to  a  linear,  exponential,  or  an  S-curved  func¬ 
tion  of  time.  A  second  form  is  decomposition,  where  a  time 
series  is  broken  into  (1)  its  long-term  trend.  (2)  its 
cyclical  component — long-term  swings  about  the  trend,  (3) 
its  seasonal  component — regular  periodic  fluctuations  in 
demand,  and  (4)  its  random  component — unpredictable  fluctua¬ 
tions.  Other  extrapolative  forecasts  are  exponential 
smoothing  (a  method  "based  on  the  weighted  sum  of  past 
observations"),  and  Box-Jenkins  (a  systematic  approach  that 
finds  the  best  particular  method  and  weighted  sums  of  a 
choice  of  models)  (11:94).  The  final  approaches  are  Bayes¬ 
ian  methods  that  permit  human  intervention  in  a  modified 
exponential  smoothing  method  to  account  for  known  changes  in 
the  environment  (11:94). 

Causal  and  Structural  Models.  The  last  group  of  fore¬ 
casting  methods  Fildes  mentions  are  the  causal  models 
(11:95).  These  methods  are  generally  more  sophisticated  and 
data  intensive. 

Fildes  identifies  a  number  of  basic  causal  forecasting 
methods:  regression,  simultaneous  systems,  simulation, 

input-output,  and  cross-impact  analysis  (11:95).  The  first 
method,  regression,  involves  finding  a  relationship  between 
demand  and  some  independent  variable  presumed  to  cause  the 
demand  (11:95).  Related  to  this  method  are  simultaneous 
system  models  which  predict  several  dependent  variables  from 
several  independent  variables  (11:95).  The  third  modeling 
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technique  is  simulation.  Simulation  models  are  more  con¬ 
cerned  with  relationships  and  less  concerned  with  strict 
mathematical  analysis  (11:95).  Another  type  of  model  is  the 
input-output  model  which  capitalizes  the  fact  that  for  a 
given  output,  a  certain  input  must  be  available  (11:95). 

Once  the  output  is  known  the  input  level  can  be  obtained. 

The  final  model  is  called  cross-impact  analysis  (11:95). 

Its  rationale  is  that  demands  are  based  on  the  occurrence  of 
certain  events.  If  the  probabilities  of  those  events  oc¬ 
curring  can  be  estimated,  then  demands  can  be  forecast. 

This  type  of  analysis  can  be  carried  out  through  several 
layers  of  events  and  probabilities  (11:95). 

Mode  1  Summary .  The  choice  of  which  method  should  be 
used  in  a  given  circumstance  depends  on  numerous  factors. 

In  fact,  "no  one  method  can  be  relied  on  to  produce  the 
'best'  forecasts  in  all  circumstances"  (11:96).  Table  1 
gives  Fi Ides'  opinion  of  the  relative  strengths  and  weak¬ 
nesses  of  the  different  methods.  Fildes  judges  each  fore¬ 
casting  method  according  to  seven  parameters:  (1)  the 
relative  number  of  data  points  required  to  use  the  method 
(DATA);  (2)  the  statistical  basis  of  the  method,  is  it 
intuitive  or  not?  (BASIS);  (3)  the  difficulty  of  setting  up 
the  method  (SET-UP);  (4)  the  relative  ease  of  using  the 
method  (EASE);  (5)  the  comprehensibility  of  the  method — how 
readily  managers  can  understand  the  intricacies  of  the 
method  (COMP.):  (6)  the  assessabi 1 i ty  of  the  method — how 
easily  the  model  can  be  evaluated,  a  measure  of  how  specific 
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the  model  is  (ASSES.);  and  (7)  the  reported  effectiveness — 
how  much  the  model  has  been  compared  to  and  evaluated 
against  other  models  in  the  literature  (EFFECT). 
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Time  Series 

Time  series  analysis  uses  past  data  to  predict  the 
future  (29:41).  By  determining  how  the  data  changed  over 
time  in  the  past,  the  future  may  be  predicted  (29:41).  Data 
behavior  may  be  broken  into  five  interacting  components.- 
" levels .  trends,  seasonal  variations,  cyclical  variations, 
and  random  variations" (  29:41).  These  are  illustrated  in 
Figure  1.  The  first  component  is  the  level  component.  This 
is  the  "central  tendency  of  a  time  series  at  any  given 
t ime " ( 29 : 41 ) .  The  trend  is  the  smooth  line  that  represents 
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Time  (years! 


the  long  term  pattern  of  the  series  (29:41).  (It  should  be 
noted  that  trend  and  level  components  are  essentially  mutu¬ 
ally  exclusive.)  Seasonal  variations  are  fluctuations  above 
and  below  the  trend  or  level  line  which  repeat  from  season 
to  season.  Normally  thought  of  as  being  annual,  some  sea¬ 
sons  are  weekly  or  even  daily.  Cyclical  variations  are 
those  that  occur  over  longer  time  periods  than  seasonal 
variations  and  tend  to  be  less  consistent  (29:41).  They 
could  be  construed  to  be  related  to  long-term  business 
cycles  if  appropriate.  Random  variations  are  those  in  the 
data  which  cannot  be  accounted  for  otherwise  and  have  no 
identifiable  pattern. 

Exponential  Smoothing 

Of  the  various  forecasting  techniques,  exponential 
smoothing  offers  one  of  the  best  combinations  of  accuracy 
and  low  data  processing  and  retention  requirements.  This 
conclusion  is  supported  in  full  by  the  literature.  Brown 
comments  that  accuracy,  simplicity  of  computation,  and 
flexibility  to  adjust  the  rate  of  response  were  three  of  the 
most  "frequently  important"  criteria  of  a  forecasting  method 
and  noted  that  exponential  smoothing  benefited  from  these 
qualities  (6:91).  Peterson  and  Silver  point  out  that  be¬ 
cause  of  the  limited  data  and  calculation  requirements  of 
simple  exponential  smoothing  models  they  are  well  suited  to 
handling  thousands  of  items  (25:111).  Tersine  notes  "the 
major  advantage  of  (simple  exponential  smoothing]  is  that 


the  effect  of  all  previous  data  is  included  in  the  previous 
forecast  figure,  so  only  one  number  needs  to  be  retained  to 
represent  the  demand  history"  (29:53).  Blazer  notes  that 
"time  series  models  are  the  most  practical  models  for  fore¬ 
casting  Air  Force  EOQ  items.”  making  special  mention  of 
single  (or  simple)  exponential  smoothing  (4:22).  In  fact. 
Patterson  notes  that  the  SB3S  uses  a  modified  exponential 
smoothing  model  (23:3) .  Since  exponential  smoothing  is  con¬ 
sidered  so  useful,  a  more  detailed  explanation  of  how  1 1. 
works  is  in  order. 

Moving  Averages.  Many  derivations  of  the  exponential 

smoothing  formula  use  moving  averages  as  a  starting  point 
(6;  18;  30).  The  bases  of  moving  averages  and  exponential 
smoothing  are  lacking  in  statistical  rigor  (30:55).  but  they 
are  still  very  useful  and  intuitively  valid  in  spite  of 
their  lack  of  a  rigorous  mathematical  foundation.  The  idea 
behind  moving  averages  is  to  forecast  based  on  an  average  of 
recent  actual  data  (30:55).  The  average  of  a  predetermined 
number  of  recent  data  points  is  used  to  forecast  the  next 
period's  demand.  The  average  "moves"  because  as  newer  data 
become  available,  the  oldest  data  points  are  eliminated  from 
consideration;  thus,  the  forecast  reflects  only  the  most 
recent  values  (30:55).  The  moving  averages  formula  is: 


Ft.  +  1  “(Xt  +  Xt-1  +  ’  •  '  +  Xt-N+l’/N 
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(2  X: >/N 
i  =  t  -  N  + 1 
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or 


-  [ (Xt )  -  1  (X  t_N+1 ) ]/N  +F.  (Ill 

(18:47-48) 

where  N  Is  the  number  of  periods  in  the  forecast.  X  is  the 
demand  in  a  period.  F  is  the  forecast  for  a  period,  t  is  the 
current  period,  and  t+1  is  the  next  period. 

It  is  apparent  from  these  formulas  that  the  greater  the 
value  of  N,  the  more  stable  the  forecast.  The  smaller  the 
value  of  N.  the  faster  the  forecast  reacts  to  changes 
(18:46).  In  other  words,  a  large  value  of  N  would  be  appro¬ 
priate  for  forecasting  stable  demand  with  relatively  large 
random  fluctuations,  because  the  forecast  would  not  follow 
the  random  variations,  but  would  smooth  the  fluctuations. 

In  contrast,  a  small  value  of  N  would  react  quickly  to 
changes  in  demand  and  follow  changes  in  demand  quickly.  In 
both  cases,  moving  averages  forecasts  would  tend  to  lag  or 
follow  the  demand  pattern  (18:46). 

According  to  Meredith,  however,  exponential  smoothing 
has  two  major  advantages  over  moving  averages:  smaller  data 
requirements  and  greater  ease  of  computation  (20:91).  Yet. 
exponential  smoothing  balances  the  need  to  smooth  the  fore¬ 
casts  with  a  number  of  periods  of  historical  data,  and  the 
need  to  forecast  new  trends  with  recent  data. 

Exponential  Smoothing.  The  basic  formula  for  exponen¬ 
tial  smoothing  can  be  developed  from  (11).  Assuming  the 
value  of  X _jsj  +  ^  is  unavailable,  a  reasonable  approximation 
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for  stationary  data  would  be  F^,  the  forecast  for  the  prev¬ 


ious  period.  Substituting  into  [11]  results  in 

Ft+1  =  Xt/N  ~  Ft/N  +  Ft 
Equation  12  can  be  rewritten  as 


[12] 

(18:49) 


F*.  .,  =  (l/N)Xt  +  ( 1-1/N)  [13] 

(18:49) 

This  results  in  a  forecast  derived  from  weighting  Xt  (the 
most  recent  data  point)  with  1/N  and  weighting  Ft  (the  most 
recent  forecast)  with  1-1/N.  Substituting  a  smoothing 
coefficient,  a.  for  1/N.  where  0  i  a  <  1,  [14]  becomes 


F  .  ...  =  (a)  Xa.  +  ( 1 -a ) Ft  [14] 

(18:49) 

Seppanen,  referenced  by  Patterson,  observes  that  in  the  SBSS 
a  =  N/(.365+N)  for  recent  demands  where  N  is  the  number  of 
days  since  the  last  semiannual  adjustment  (23:3).  Exponen¬ 
tial  smoothing's  reduced  data  storage  requirements  can  be 
readily  seen:  only  the  most  recent  demand,  the  most  recent 
forecast  and  a  must  be  retained.  Rewriting  [14]  yields 


or 


t  +  1 


Ft  +  a(Xt  -  Ft) 


[15] 

(18:501 


Ft  +  1  Ft  +  aet 


[16] 

(18:50) 


where  efc  is  the  forecast  error  for  period  t.  So.  the  new 
forecast  for  period  t  +  1  is  actually  the  previous  period's 
forecast  plus  a  fraction  of  the  error  between  forecast  and 
demand.  Intuitively,  this  seems  appropriate  that  a  good 
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estimator  of  the  future  would  be  the  last  forecast  plus  a 
portion  of  the  error  in  making  that  forecast  (18:48-50). 

Exponential  Smoothing  with  Seasonal  Factor . 

Winters  (one  of  the  original  developers  of  exponential 
smoothing)  expanded  the  usefulness  of  the  model  by  incor¬ 
porating  a  seasonal  demand  factor.  He  started  by  assuming 
that  any  seasonality  was  proportional  to  the  demand  of  that 
period  rather  than  being  additive,  or  independent  of  demand 
(32:327).  The  first  equation  of  the  Winters'  model  can  be 
expressed  as : 


St  =  (a)Xt/It_L  +  (l-a)Ffc_1 


[171 

(32:3281 


where  L  is  the  number  of  periods  in  one  seasonal  cycle,  (for 
example,  for  a  yearly  cycle  L  would  be  12)  and  St  is  the 
"estimate  of  the  expected  deseasona 1 ized  sales  rate  in 
period  t“  (32:328).  The  estimated  seasonal  factor  for  the 
period  t.  It  .  can  be  calculated  using 


It  =  (0>Xt/St  ♦  d-B)It_L 


[18] 

(32:328) 


where  0  is  the  seasonal  smoothing  factor,  which  must  be 
chosen  empirically.  The  seasonal  factor.  1^.  .  is  chosen  so 
that  like  periods  are  forecast:  the  last  corresponding 
period  must  be  used  in  the  calculations.  For  example,  one 
must  use  last  year's  June  I  ^  value  to  make  this  year's  June 
forecast  (32:328).  The  actual  forecast  is  made  using: 


Ft+M  “  St  * rt-L+M 
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(32:328) 


where  M  is  the  number  of  periods  into  the  future  to  fore¬ 
cast  (32:327-329). 

Exponential  Smoothing  with  Trend  and  Seasonal 
Factors .  The  previous  model  does  not,  however,  compensate 
for  any  type  of  trend.  Winters  also  developed  a  trend  and 
seasonal  model  where: 


=  (a)Xt/It. 


+  ( 1-a) (F 


t-1 


bt-}J 


32:329 


[20] 

18:72) 


The  b^_^  factor  accounts  for  an  additive  trend.  The  season¬ 
al  factor  remains  unchanged  from  [18]: 


I.  -  ( 8 )  Xt/St  +  (H)L.  [21] 

(32:330;  18:72) 

The  trend  factor  resembles  [19]  and  [20]  in  its  structure: 
(32 : 330) 


bf  -  (D(S.-St  ,)  +  (1-Db,  ,  [22] 

(32:330:  18:72) 

The  smoothing  constant.  T.  must  be  chosen  arbitrarily.  The 
final  forecast  becomes 


Ft+M  ^St  +  bt*M)*Xt-L+M 


[23] 

(32:330:  18:73) 


To  obtain  a  forecast  using  either  of  these  methods.  St 
is  first  calculated  using  because  1^.  cannot  be  deter¬ 

mined  until  St  is  calculated.  Then  It  and  bt,  if  required 
can  be  obtained.  Finally  Ft+j^  can  be  found  from  [19]  or 
[23]  (32:329-330:  18:72-73). 

Useful  as  Winters'  method  is,  it  has  one  major  draw¬ 
back:  "it  requires  three  smoothing  parameters"  (18:79). 
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This  greatly  increases  the  data  requirements  and  computer 

time  required  to  generate  forecasts.  In  response  to  this 

concern.  Zehna  and  Taylor  noted  as  long  ago  as  1975  that 

computer  processing  times  should  not  “be  a  significant 

factor  in  the  selection  of  forecasting  techniques"  because 

of  the  great  increases  in  computer  power  and  speed  (33:38) . 

Ammer  adds  that  in  general: 

Exponential  smoothing  is  not  a  perfect  technique 
for  forecasting  demand.  Like  other  statistical 
techniques,  it  is  based  on  past  demand  and  obvi¬ 
ously  cannot  allow  for  new  and  unpredicted  changes 
in  demand.  In  addition,  it  may  help  increase  in¬ 
ventories  because  .  .  .  small  shifts  in  weekly 

demand  can  be  magnified  into  substantial  shifts  in 
order  points  when  trend  is  taken  into  account. 

Despite  its  weaknesses,  however,  exponential 
smoothing  is  an  excellent  technique,  especially 
when  used  in  conjunction  with  electronic  data  pro¬ 
cessing  to  control  order  points  of  thousands  of 
low-value  C  items.  (2:305) 


Box-Jenkins 

A  more  sophisticated  methodology  for  forecasting  which 
may  outperform  exponential  smoothing  was  developed  by  Box 
and  Jenkins  to  create  a  class  of  forecast  models  known  as 
ARIMA  models.  This  method  is  based  on  two  major  premises 
(30:171).  First,  the  data  are  assumed  to  be  in  the  form  of 
a  time  series,  in  which  observations  are  taken  over  a  period 
of  time,  and  the  forecast  is  needed  for  some  future  period 
in  time  (30:171-172).  Second,  the  data  are  assumed  to  be 
following  a  pattern  with  an  obscuring  random  component 
(30:171-172).  The  objective  of  the  forecaster,  then,  is  to 
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identify  the  pattern  and  use  it  to  make  forecasts  (30:171- 
172)  . 

Box-Jenkins  can  be  considered  a  process.  The  follow¬ 
ing  iterative  procedure,  diagrammed  in  Figure  2.  is  recom¬ 
mended  by  Box  and  Jenkins  to  develop  a  forecast  (5:19).  The 
essential  idea  is  that  an  initial  class  of  model  is  selected 
based  on  certain  factors  (to  be  discussed  in  greater  detail 
later),  and  a  model  is  tentatively  identified  (5:18).  The 
model  parameters  are  then  estimated  and  the  model  is  eval¬ 
uated  against  certain  diagnostic  checks  (5:19).  Based  on 
the  evaluation,  either  the  model  is  used,  or  another  model 
is  tested  (5:19).  Eventually,  a  model  meeting  the  diagnos¬ 
tic  requirements  is  developed  and  the  model  is  selected  to 
make  forecasts  (5:19). 


Figure  2.  Iterative  Model  Building  Process  (5:19) 
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The  Box-Jenkins  methodology  incorporates  two  major 
types  of  models  which  together  make  up  the  class  of  autore¬ 
gressive  integrated  moving  average  processes  (ARIMA 
processes)  (5:8).  In  order  to  better  understand  the  ter — 
mi no  logy  involved,  the  backshift  operator  is"  first  defined 
as 


Bz 


t 


ry 

n-i 


The  more  general  case  is 


[24] 

(5:8) 


Bmz* 


't-m 


[25] 

(5:8) 


Another  important  operator  is  the  backward  difference 
operator.  6.  which  may  be  defined  by  the  following  equation.- 


6zt  =  (1-B)zt  -  zt  -  zt_! 


[26] 

(5:8) 


These  operators  are  computational  artifacts  which  make  the 
mathematical  analysis  more  straightforward  and  permit  the 
simplification  of  a  number  of  equations. 

Autoregressive  Models.  A  very  useful  model  in  the  Box- 
Jenkins  methodology  with  practical  applications  is  the 
autoregressive  model.  In  this  model  the  "current  value  of 
the  process"  is  a  function  of  the  previous  values  of  the 
process  and  some  additional  factor  or  shock.  (5:9).  If 
the  values  of  the  process  at  times  t,  t-1 ,  t-2  .  .  .  are 

represented  by  zt .  zt-2  ‘  •  ■  zt.  2t-l'  zt-2  can 

be  the  differences  of  each  zfc  from  the  mean  u.  as  in  zfc  =  zt 
-  u .  Then 
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“  £  0  l^t  — 1  ^2zt-2  + 


■  +  0pZt-p  +  at 


[27] 

(5:9) 


is  an  autoregressive  ( AR )  process  of  order  p.  Using  the 
autoregressive  operator  of  order  p 


0B  =  1  -  f^B  -  <zS2B2 


a  bp 

p 


[28] 

(5:9) 


the  AR  model  may  be  written 


d(B)zt  =  afc 
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(5:9, 


zt  -  0(B) afc 


[30] 

(5:9) 


where 


0(B)  -  <zi  i(B) 


[31] 
(5 : 10) 


Moving  Average  Models.  A  second  model,  the  moving 
average  process,  is  based  on  the  at's,  the  error  terms.  A 
moving  average  (MA)  process  of  order  q  may  be  written 


zt  31  at  ®i®t  —  1  ®2at  —  2 


6,^4.  „ 

q  t-q 


(5:10) 


If  a  moving  average  operator  is  defined  as 


9(B)  -  1  -  ©1B  -  92B2 


vq 


[33: 

(5:10 


then  the  MA  model  simplifies  to 


zt  =  0(B)at 


1341 

(5:10) 
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Zt(l  -  0B) 


-1 


135] 

(5:10) 


or 

at  -  ztir 


[361 

(5:50) 


where 


ir(B)  =  ( 1 -6B ) 


[37] 

(5:50) 


Autoregressive  -  Moving  Average  Models.  If  the  data 
support  it.  a  combination  of  AR  and  MA  terms  is  necessary. 
The  ARMA  model  is  simply  an  AR  and  an  MA  model  combined.  It 
may  be  written  as  follows: 


zt  *  c>lzt-l  + 


+  *pzt-p  +  at  " 

9lat-l  '  •  •  ■  "  6qat-q 


or  in  shorthand 


[381 

(5:11) 


d(B)zt  -  0(B) at  [391 

(5:11) 

Although  the  AR.  MA ,  or  combined  models  suit  stationary  time 
series  (those  that  vary  about  a  constant  mean),  frequently 
the  data  call  for  a  non-stat ionary  model  that  will  account 
for  changes  in  the  mean  of  the  process.  By  introducing  the 
"generalized  autoregressive  operator"  ?(B)  and  the  station¬ 
ary  operator.  e)(B).  for  the  AR  process,  one  can  write  (5:11) 

<P(  B)zt  -  <z!  (B)  ( 1-B)dzt  =  9  ( B)  a.  [40] 

(5:11) 

where  the  term  (1-B)  accomplishes  the . d i f f erenc i ng .  In 
other  words,  an  ARMA  process  is  differenced  by  creating  a 
new  series  with  the  differences  between  the  value  of  one 
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period  and  the  value  of  the  previous  period.  The  new  series 
is  a  more  generalized  model.  This  gives  a  general  model 


#  ( B )  w^.  =  © (B) at  [41] 

(5:11) 

where 

wt  =  6d zt  [42] 

(5:11) 


Thus,  by  taking  d  differences  .of  the  process  to  make  it 
stationary,  nonstationary  behavior  may  be  modeled.  This 
last  process  described  is  called  an  autoregressive 
integrated  moving  average  (ARIMA)  process  of  order  (p.  d.  q) 


(5:11). 


Figure  3  illustrates  how  the  ARIMA  process  "filters" 
random  errors  through  the  three  operations  to  yield  the  time 
series . 

Seasona 1 itv .  One  major  characteristic  of  a  seasonal 
time  series  is  that  every  s  (for  annual  seasonality  of 
monthly  data,  s  -  12)  periods  a  similar  observation  may  be 
made.  Another  characterist ic  is  that  successive  months  (in 
the  case  of  annual  seasonality)  will  also  be  similar.  For 
example,  an  observation  made  in  October  will  be  similar  to 
one  made  the  previous  October  and  all  other  October  observa¬ 
tions.  Additionally,  the  October  observation  will  be 
similar  to  September's  and  November's  observations  for  the 
same  year.  (5:303-304) 


29 


0(B) 


Moving 
Average 
Fi Iter 


fiJ  1  CB) 


Stationary 
Auto¬ 
regressive 
Fi 1  ter 


St 

Time  Series 


Non¬ 
stationary 
Summation 
Fi  Iter 


Random  Error 


Figure  3.  Filtering  Process  to  Convert  White  Noise  to 
Time  Series  (5:12) 


A  model  of  the  form 


*(B°)6^zt  =  6(B  )at 


(431 
(5 : 304) 


may  be  written  where  s  -  12,  6S  **  1  -  Bs ,  and  4(BS)  and  0 
(Bs)  "are  polynomials  in  Bs  of  degrees  P  and  Q,  respective¬ 
ly.  and  satisfying  stationarity  and  invert ibi 1 ity  condi¬ 
tions"  (5:304).  A  stationary  model  randomly  varies  about 
one  mean  (5:7)  and  has  the  condition  that^(B)  is  within  the 
range  -1  to  +1  (5:51).  An  invertible  model  is  slightly  more 
difficult  to  explain.  Because  the  moving  average  model  de¬ 
pends  on  the  error  terms  for  its  value  of  .  if  the  ab¬ 
solute  values  of  the  individual  0’s  are  greater  than  one. 
the  weights  of  the  at ' s  will  increase.  leading  to  a 
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nonstationary  condition  (5:50).  Thus  the  condition  of 
invert ibi  1  ity  is  imposed,  that  the  sum  of  the  01BJ  =  ir( B)  is 
within  the  range  -1  to  +1  (5:51). 

Having  defined  and  described  Box-Jenkins  models  the 
methodology  behind  applying  them  will  be  discussed  in 
Chapter  3. 

Previous  Studies 

The  management  of  seasonal  items  has  received  much 
attention  in  the  literature  (6;  21;  25:  29:  30).  although 
few  military  studies  have  been  done.  One  of  these  was  the 
Air  Force  Institute  of  Technology  thesis  by  Gil  loth,  Ohl , 
and  Wells  in  1979  titled  "An  Evaluation  of  Seasonality  in 
the  United  States  Air  Force  Medical  Material  Management  . 
System"  (MMMS)  (14).  After  examining  the  demand  patterns  of 
1886  items  (14:36).  they  determined  that  seasonality  was 
evident  in  25  -  35%  of  the  medical  supplies  examined,  while 
only  about  8  -  9%  of  the  items  were  seasonal  over  both  the 
years  tested.  They  concluded  that  double  exponential 
smoothing  was  the  best  forecasting  method  tested,  using  mean 
squared  error  as  the  selection  criteria  (14:28.58).  Other 
forecasting  methods  examined  were  the  current  MMMS  system 
(12  month  moving  average),  and  adaptive  response  rate  ex¬ 
ponential  smoothing  (a  method  in  -which  the  smoothing  con¬ 
stant  is  automatically  varied  as  the  demand  pattern  changes) 
(14:43).  One  of  their  recommendations  was  for  increased 
data  retention  to  benefit  future  studies,  since  they  only 
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had  two  years  of  non-consecutive  data  with  which  to  work 
(14:59).  This  lack  of  data  limited  considerably  the  gen¬ 
era  1  izabi  1  ity  of  their  research  (14:20). 

Other  studies  include  one  done  by  Orchowsky  in  1985  on 
clothing  items,  in  which  some  items  investigated  showed 
seasonal  demands,  but  most  items  did  not  (22:30-31).  This 
study  examined  4563  individual  national  stock  numbers, 
managed  by  the  Defense  Logistics  Agency,  which  were  grouped 
into  300  Procurement  Grouping  Codes  (PGC)  (22:3).  (A  PGC  is 
a  group  of  like  clothing  items  differentiated  only  by  size 
(22:3)).  Orchowsky  used  the  autocorrelation  function  (ACF) 
for  12  lags  to  determine  if  a  random  or  seasonal  pattern 
existed  (22:10).  Briefly  "the  ACF  is  a  measure  of  the 
relationship  (correlation)  of  the  time  series  with  itself, 
lagged  by  some  number  of  time  periods"  (22:10).  For  example 
data  points,  taken  monthly.  12  periods  apart,  with  a  high 
autocorrelation  would  indicate  a  strong  correspondence  on  a 
seasonal  basis — seasonal  data. 

Fischer  and  Gibson,  in  1971.  examined  single,  double 
(with  a  trend)  and  triple  (with  a  quadratic  component) 
exponential  smoothing  and  a  12  month  moving  average  (12:9). 
Their  sample  consisted  of  34  EOQ  items  from  Wr ight-Patterson 
AFB  Base  Supply.  Fischer  and  Gibson  were  disappointed  in 
the  performance  of  triple  exponential  smoothing,  expecting 
it  to  perform  better  than  it  did  (12:67).  They  concluded 
that  an  inadequate  database,  with  only  22  months  of  data 
available,  prevented  the  determination  of  acceptable 
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smoothing  constants  (12:71).  The  smoothing  constants  were 
"optimized"  by  testing  several  values  for  the  minimum  sum  of 
squared  errors  (12:45).  Statistically,  the  four  models  had 
the  same  mean  square  error,  with  or  without  "optimized" 
smoothing  constants  (12:49). 

A  fourth  study,  undertaken  by  Bittel  and  Gartner  in 
1982,  examined  demand  for  consumable  items  at  the  depot 
level  (3).  Their  analysis  of  800  line  items  (3:42)  used  the 
following  models  (3:56): 

1.  Naive  (the  next  period  forecast  is  the  last 

period's  actual  value) 

2.  Simple  moving  average  (4.  8,  and  12  periods) 

3.  Double  moving  average  (4.  8.  and  12  periods) 

4.  Single  exponential  smoothing  (a=0.2  and  0.8) 

5.  Focus  forecasting  (a  "multi-model  technique  which 
employs  simplistic  forecasting  assumptions  and 
computer  simulation  to  forecast  demand"  (3:49- 
50)  ) 

6.  Simple  regression  (a  method  of  fitting  the  line 
with  the  least  mean  squared  error  to  the  data) 

7.  S-curve  analysis  (a  model  used  to  describe  and 
"estimate  the  life  cycle  of  technologies  and  prod¬ 
ucts.  An  S-curve  has  a  slow  start,  a  rather  steep 
growth,  and  a  saturation  that  comes  after  some 
period  of  time"  (19:169)) 

8.  Exponential  growth  (a  model  which  describes  con¬ 
stant  growth  rates  over  long  periods  of  time 
(19:171) ) 

9.  Eclectic  forecasting  (similar  to  focus  forecast¬ 
ing.  but  uses  more  complex  techniques  (3:50)) 

Although  Bittel  and  Gartner  were  looking  at  a  different 

echelon  of  demand,  they  noted  that  simple  exponential 

smoothing  gave  the  lowest  mean  average  deviation  and 
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variance  of  all  the  methods  tested.  However,  they  observed 
that  the  method  was  statistically  equivalent  to  single  [sic] 
moving  averages  with  four  and  eight  month  periods  (3:79). 
Additionally,  they  noted  that  the  data  "showed  some  normal 
demand  patterns  with  considerable  random  variation"  (3:75). 

Another  study  which  indirectly  examined  seasonality  was 
one  by  Gertcher  in  1982  (13).  As  part  of  his  research  into 
different  methods  of  forecasting  demand  for  consumable 
items,  he  used  triple  exponential  smoothing,  which  assumes  a 
quadratic  shape  to  a  plot  of  the  historical  demand  data 
(13:20).  His  major  difficulties  in  using  triple  exponential 
smoothing  were  poor  initialisation  constants  and  the  lack  of 
fit  of  the  demand  data  to  the  assumed  second  order 
polynomial  (13:44).  In  other  words,  the  data  did  not  fully 
exhibit  the  quadratic  pattern  of  the  model  he  was  using. 

The  consensus  of  most  of  these  studies  is  summed  up  in 
a  study  on  demand  variance  by  Blazer,  who  noted  that  the 
nemesis  of  the  analysis  of  competing  forecasting  methods  is 
a  lack  of  sufficient  demand  data  (4:3).  He  notes  that  10 
years  of  data  is  desirable.  Yet  Gil  loth,  Ohl .  and  Wells 
used  data  from  calendar  years  1975  and  1977  (14:20); 
Orchowsky  was  able  to  obtain  33  months  of  data  (22:2); 

Zehna  and  Taylor  had  only  8  quarters  of  data  (33:22).  Each 
of  these  efforts  suffered  from  the  lack  of  data  maintained 
by  the  SBSS . 
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Determination  of  Seasonality 


Peckham  gives  three  criteria  for  manually  determining 
if  seasonal  trends  exist: 

(1)  The  peak  demand  should  be  substantially  higher 
than  the  random  fluctuations  or  "noise"  in  the 
demand . 

(2)  The  peak  demand  must  occur  during  the  same  time 
period  each  year. 

(3)  The  reason  for  the  peak  must  be  known.  (25:40-41) 
Based  on  this,  a  plot  of  the  data  would  frequently  give  an 
indication  of  whether  or  not  an  item  exhibited  seasonal 
demands.  A  strong  indication  of  seasonality  would  be  a 
pattern  of  demands  generally  repeating  itself  over  several 
years,  with  peaks  and  valleys  being  considered  significant 
if  greater  than  two  standard  deviations  from  the  mean  de¬ 
mand.  This  would  give  the  demand  spikes  a  4 %  chance  of 
being  considered  significant  when  they  actually  were  not, 
assuming  a  normal  distribution  of  demands. 

However,  manually  determining  seasonality  in  this  fash¬ 
ion  does  not  lend  itself  to  a  statistically  verifiable 
analysis.  Without  a  definition  of  "substantially  higher" 
different  analysts  may  arrive  at  different  conclusions  re¬ 
garding  the  same  data.  A  method  relying  less  on  intuition, 
the  ACF  as  used  by  Orchowsky  (23).  could  be  readily  used  to 
detect  seasonality  in  data. 

To  better  understand  autocorrelation,  simple  corre¬ 
lation  will  be  defined  first,  as  in  a  discussion  by  Wheel¬ 
wright  and  Makridakis  (30:116).  Simple  correlation  is  a 
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measure  of  the  relationship  between  two  variables.  A 
correlation  coefficient,  r,  which  ranges  from  -1  to  +1.  is 
used  to  measure  the  strength  of  the  relationship.  With 
perfect  correlation,  r  would  be  either  ±1:  with  no  correla¬ 
tion  r  would  be  0 .  A  value  of  +1  is  assigned  for  a  positive 
correlation,  while  a  value  of  -1  is  assigned  for  a  negative 
correlation.  In  other  words,  r  would  be  +1  between  two 
identical  sets  of  numbers,  and  -1  between  a  set  of  negative 
numbers  and  the  absolute  values  of  those  same  numbers.  In 
addition,  an  r  of  0  implies  there  is  no  relationship  between 
the  two  variables,  that  is:  they  are  completely  independent 
and  a  change  m  one  has  absolutely  no  effect  on  the  other. 
(30:116.  174) 

An  autocorrelation  coefficient  extends  this  same  idea 
to  "values  of  the  same  variable  but  at  different  time 
periods"  (30:174).  In  this  case,  the  data  under  considera¬ 
tion  are  made  into  several  different  data  sets,  each  treated 
as  a  different  variable.  The  first  data  set  consists  of  the 
original  data.  The  second  data  set  is  the  same  data,  with 
the  first  value  removed,  and  all  others  moved  up  one  spot. 
Thus,  the  original  value  for  period  1  is  dropped,  the  origi¬ 
nal  value  for  period  2  becomes  the  new  value  for  period  1. 
the  original  value  for  period  3  becomes  the  new  value  for 
period  2  and  so  on.  Then  this  new.  second,  data  set  is  used 
to  generate  a  third  data  set  using  the  same  principles.  The 
second  data  set's  value  for  period  1  is  dropped,  the  second 
data  set's  value  for  period  2  becomes  the  newest  value  for 
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period  1.  the  second  data  set's  value  for  period  3  becomes 
the  newest  value  for  period  2  and  so  on.  Thig  is  done  as 
long  as  data  is  available  to  continue  creating  new  data 
sets.  ( 30 : 174-175 ) 

Wheelwright  and  Makridakis  continue  by  observing  that 
each  data  set  can  now  be  considered  a  new  variable,  and  a 
correlation  coefficient  can  be  calculated  between  each  of 
the  data  sets  (30:175).  The  correlation  coefficient  tells 
how  data  a  certain  number  of  time  periods  apart  is  cor¬ 
related.  The  correlation  between  the  original  data  set  and 
the  second  data  set  shows  how  data  one  period  (lag)  apart 
are  correlated  and  whether  they  tend  to  move  in  the  same 
direction  or  not.  The  same  can  also  be  done  with  the  rest 
of  the  data  sets  (30:175).  The  correlation  between  the 
original  data  and  the  lagged  data  is  obtained.  Since  the 
correlation  is  between  the  data  and  itself,  one  or  more 
periods  removed,  the  correlations  are  called  autocorrela¬ 
tions  (30:175-176).  The  ACF  is  useful  in  identifying  sea¬ 
sonal  data  because  it  can  be  used  to  identify  at  which  lags 
the  correlations  are  largest.  If  the  ACF  for  the  12th  lag 
for  monthly  data  (data  points  12  months  apart)  is  statisti¬ 
cally  significant  then  the  data  can  be  reasonably  presumed 
seasonal  (30:176). 

Forecast  Performance  Measures 

A  number  of  methods  are  available  to  determine  how 
"good"  a  particular  forecasting  method  is  at  forecasting. 
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The  basic  assumption  of  any  forecasting  technique  is  that 
there  is  a  fundamental  pattern  to  the  observations,  plus 
some  random  error  or  fluctuations.  The  object  of  fore¬ 
casting  then,  is  to  minimize  those  errors.  Error  is  usually 
defined  as  the  difference  between  the  actual  value  and  the 
forecasted  value. 

A  method  of  evaluating  those  errors  is  necessary.  A 
simplistic  method  of  analyzing  error  would  be  to  average  all 
the  errors,  but  since  positive  and  negative  errors  may 
cancel  each  other  out.  this  could  lead  to  an  erroneous 
assumption  that  the  overall  errors  are  small.  An  improve¬ 
ment  would  be  to  average  the  absolute  value  of  the  errors, 
but  this  would  weight  all  errors  the  same.  Yet  another 
alternative  would  be  to  average  the  squared  errors.  This 
more  heavily  penalizes  large  errors  in  favor  of  small  ones . 
An  error  of  4  would  be  counted  16  times  as  much  as  an  error 
of  1.  This  method,  the  Mean  Squared  Error  (MSE) .  will  be 
used  in  this  study. 

Conclusions 

The  SBS5  has  a  sophisticated  demand  forecasting  system 
designed  to  compute  an  average  daily  demand  rate,  which  is 
used  to  establish  stockage  levels.  The  system  can  be 
thought  of  as  a  modified  exponential  smoothing  model.  An 
examination  of  a  variety  of  forecasting  methods  revealed  a 
number  of  potential  forecasting  methods.  However,  experts 
tend  to  agree  that  if  many  items'  demand  must  be  forecast. 
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exponential  smoothing  has  several  advantages  over  other 
methods,  especially  its  minimal  computational  and  data 
storage  requirements.  Exponential  smoothing  is  a  very 
simple  method  of  forecasting,  although  the  addition  of  terms 
to  compensate  for  trend  and  seasonality  added  considerable 
complexity  to  the  model.  In  spite  of  this,  exponential 
smoothing  remains  the  model  of  choice  for  forecasting  large 
numbers  of  items  because  of  its  modest  requirements.  The 
Eox-Jenkins  methodology  is  an  powerful  wav  of  generating  the 
best  time  series  or  ARIMA  model  to  fit  a  given  data  set. 
Although  relatively  difficult  to  use.  it  offers  confidence 
intervals  for  both  parameters  and  estimates,  in  contrast 
with  the  intuitive  exponential  smoothing  techniques.  In 
addition,  the  Box-Jenkins  methodology  is  quite  data  inten¬ 
sive. 

Little  military  research  has  been  performed  in  examin¬ 
ing  how  well  seasonal  forecasting  methods  perform  with 
consumable  SBSS  items.  Previous  studies  that  did  attempt 
this  had  mixed  success,  primarily  because  of  the  lack  of 
sufficient  data. 

The  next  chapter  discusses  how.  the  actual  research  was 
conducted.  It  also  describes  how  the  data  were  obtained  and 
ana lyzed . 
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Ill .  Methodology 


Introduction 

This  chapter  discusses  the  data  collected  for  the  re¬ 
search  and  how  it  was  obtained,  manipulated,  and  finally 
used.  The  chapter  also  mentions  the  limitations  and  assump¬ 
tions  involved  in  the  study.  The  actual  research  plan  is 
highlighted  along  with  details  of  the  SBSS.  Exponential 
Smoothing  and  Box-Jenkins  methodologies. 

Data  Co  1  lection 

Population  and  Sample.  The  population  of  interest  is 
the  entire  set  of  consumable  items  used  by  the  United  States 
Air  Force.  The  sample  will  consist  of  twelve  out  of  77 
consumable  items  identified  as  seasonal  by  Langley  AFB .  VA . 
which  are  listed  in  Attachment  1.  Only  twelve  items  were 
selected  for  study  because  of  the  difficulty  in  extracting 
demand  data  from  transaction  history  tapes. 

The  data  are  not  genera  1 1 zab 1 e  to  the  entire  popul¬ 
ation  of  Air  Force  bases  because  the  selection  of  Langley  as 
a  study  base  was  non-random.  Demand  patterns  depend  on  a 
number  of  factors  including  how  and  when  certain  items  are 
ordered  at  a  facility,  the  weapons  systems  at  the  facility, 
the  mission  of  the  facility,  and  potentially  its  geographic 
1 ocat i on . 

Assumpt ions .  The  following  assumptions  were  made  re¬ 
garding  this  investigation: 
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1.  Responses  from  the  Stock  Control  Section  at 
Langley  AFB  accurately  reflect  its  situation. 

2.  Data  transmitted  to  and  from  the  Air  Force  Logis¬ 
tics  Management  Center  (AFLMC)  is  accurate  and 
comp  1 ete . 

3.  Changes  in  the  climatic  environment  of  Langley  AFB 
over  the  data  history  are  such  that  they  do  not 
significantly  impact  the  seasonality  of  the  data. 

Limitat ions .  The  following  limitations  are  apparent  in 
this  investigation: 

1.  The  limited  number  of  both  bases  investigated  and 
years  of  data  do  not  permit  genera  1 izabi 1 ity  of 
the  results  Air  Force-wide. 

2.  Variations  in  the  methods  Langley  AFB's  Stock 
Control  Section  uses  to  determine  which  items  are 
seasonal  or  not  may  make  those  items  they  deter¬ 
mine  to  be  seasonal  unrepresentative  of  other 
bases  . 

3.  No  aircraft  parts  were  identified  as  seasonal: 
only  items  such  as  cold  weather  gear  and  deicers 
were  identified  as  seasonal. 

Data  Collection.  The  data  came  from  two  sources. 

First,  Mr  Almond,  the  Stock  Control  Officer  at  Langley  AFB 
contributed  a  list  of  consumable  items  that  the  base 
considers  seasonal.  The  purpose  in  obtaining  this  list  was 
to  determine  if  selected  items  on  the  list  indeed  hold 


41 


seasonal  characteristics  in  their  demand  patterns  and  to 
develop  forecast  models  for  selected  items  from  the  sample. 

The  second  source  of  data  was  the  database  at  AFLMC . 

This  database  was  started  to  provide  historical  data  beyond 
that  normally  held  by  the  SBSS.  The  necessary  data  were  ex¬ 
tracted  from  two  types  of  records:  item  records  and  trans¬ 
action  histories.  Item  records  are  the  maior  record  on  each 
item  that,  together  with  detail  records,  collects  the  data 
required  "to  manage  most  items  under  nearly  all  circum¬ 
stances"  (8:4-64).  Item  records  were  obtained  that  reflect 
the  end  of  March  1989.  Transaction  histones  are  a  con¬ 
solidated  record  of  all  transactions  that  affect  the  inven¬ 
tory  of  an  it^m.  Unfortunately  transaction  histories  were 
only  available  for  the  dates  marked  by  an  X  in  Figure  4. 

The  extracted  data  included  national  stock  number:  nomen¬ 
clature;  expendab l 1 i ty :  repa i rabi 1 l ty :  recoverabi 1 i ty  categ¬ 
ory  code  ( ERRC ) :  cumulative  recurring  demands;  dates  of 
first  demands;  and  demands.  The  first  three  data  sets  are 
available  from  the  item  record . 
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Figure  4.  Data  Available  (X)  from  AFLMC 
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The  demand  data  were  obtained  through  analysis  of  the 
transaction  histories  for  each  item  for  the  period  of 
interest.  A  demand  was  considered  to  have  occurred  for 
every  Document  Identifier  Code  (DIC)  of  ISU  (for  issue)  and 
DUO  (for  due-out)  for  the  item.  As  a  check  of  the  ISU  and 
DOR  codes.  Type  Transaction  Phrase  and  Codes  (TTPC)  were 
used  (26).  A  TTPC  code  of  "1A"  which  decreases  the  avail¬ 
able  quantity  recorded  under  the  item  record  detail,  a  code 
of  "2D"  which  adds  and  increases  the  due-out  detail,  a  code 
of  "3Q"  which  indicates  a  post-post  (computer  off-line) 
issue  were  the  codes  considered  to  substantiate  a  request 
for  an  issue  from  stock  (8:3-134-137). 

Both  mainframe  and  personal  computers  were  used  to 
extract  the  data.  After  the  data  were  shipped  from  AFLMC  on 
magnetic  tape,  they  were  loaded  on  the  Air  Force  Institute 
of  Technology's  VAX  mainframe  computer.  One  of  several  FOR¬ 
TRAN  programs  were  then  used  (an  example  is  in  Attachment  2) 
to  extract  only  those  stock  numbers  selected  for  further 
study  and  load  the  needed  portions  of  the  transaction  re¬ 
cords  for  each  stock  number  into  a  separate  file.  Each  file 
was  then  transferred  via  modem  to  a  floppy  diskette  where  it 
was  accessible  to  an  MS-DOS  personal  computer.  Once  on  a 
diskette,  each  file  was  imported  into  Quattro.  a  spreadsheet 
program.  Then  the  database  management  function  of  Quattro 
was  used  to  extract  only  those  records  with  either  an  ISU  or 
DUO  DIC  code  and  with  a  transaction  date  within  a  given 
month.  TTPC  codes  were  manually  verified  along  with  the 
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requested  stock  number  to  ensure  the  item  requested  was  the 
item  ordered.  The  number  of  units  requested  for  each  demand 
transaction  were  summed  to  give  a  final  number  of  units  de¬ 
manded  for  each  month  for  each  stock  number.  This  figure 
was  hand-copied  and  subsequent ly‘ veri f ied .  The  resulting 
demands  for  the  twelve  stock  numbers  for  each  available 
month  are  listed  in  Attachment  3. 

One  problem  with  the  data  was  the  presence  of  a  gap 
between  December  1983  and  August  1986.  According  to  Mr. 
Miller  of  the  AFLMC  all  of  the  data  in  the  AFLMC  database 
has  a  similar  gap  (varying  in  length  at  different  bases)  be¬ 
cause  of  the  installation  of  new  base  supply  computer  sys¬ 
tems  during  that  time  period  (21).  He  stated  that  the  AFLMC 
has  "given  you  everything  we've  got"  (21).  Because  the  Box- 
Jenkins  methodology  requires  so  much  data,  it  would  be 
impossible  to  develop  any  models  on  the  relatively  con¬ 
tinuous  August  1986  through  April  1989  data. 

To  determine  the  feasibility  of  using  data  from  before 
and  after  the  gap  the  following  procedure  was  used.  The 
demand  data  from  1982  and  1983  for  each  stock  number  was 
compared  to  the  same  stock  number's  demand  data  from  1987. 
1988,  and  1989  to  statistically  determine  if  the  two  groups 
of  data  came  from  the  same  distribution.  A  Wilcoxon  Signed 
Rank  Test  was  used.  This  technique  tests  if  two  samples 
came  from  the  same  population  distribution  or  not.  The 
results  of  the  test  are  in  Table  TI.  They  indicate  that 
only  seven  of  the  twelve  items  had  demand  data  useable  for 
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Table  II.  Wilcoxon  Signed  Rank  Test  Results  for  Similar 
Distributions 


NSN 


1 

2 

3 

4 

5 

6 
7 
3 

9 

10 
11 
12 


the  Box-Jenkins  methodology  with  a  90%  confidence  level .  If 
the  P-value  was  greater  than  0.1  then  the  null  hypothesis 
that  the  two  distributions  were  the  same  could  not  be  re¬ 
jected.  For  these  seven  items  the  1982  and  1983  demand  data 
were  treated  as  if  they  came  from  the  years  1985  and  1986. 
and  combined  with  the  1987  and  1988  data.  The  1989  demand 
data  were  used  as  a  holdout  sample  for  comparison  purposes. 
Since  the  exponential  smoothing  techniques  require  less 
data,  they  were  run  using  longest  continuous  string  of  data 
available:  August  1987  through  December  1988,  with  the  1989 

demand  data  saved  as  a  hold-out  sample.  Thus,  the  exponen¬ 
tial  smoothing  forecasts  are  based  on  a  slightly  different 
database  than  the  ARIMA  forecasts. 
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Unfortunately,  as  was  noted  in  Figure  4,  the  data  was 
incomplete,  and  several  months'  data  was  lacking  in  addition 
to  the  1983-1986  gap.  This  posed  a  severe  problem  for  time 
series  analysis  since  continuous  data  is  essential  to  the 
methodology.  Peterson  and  Silver  note  that  "data  will  not 
always  be  available,  [so  one]  must  resort  to  asking  knowl¬ 
edgeable  persons  to  provide  .  .  .  guesstimates  for  the 

values  of  level,  trend,  and  seasonal  factors"  (25:120). 

They  additionally  note  that  the  forecasts  obtained  from 
incomplete  data  augmented  by  expert  opinion  must  be  care¬ 
fully  monitored  and  the  forecasting  model  updated  as  new 
data  becomes  available. 

To  fill  in  the  missing  data  points,  a  Delphi  technique 
was  used.  This  technique  involves  using  experts  to  forecast 
when  data  is  otherwise  unavailable.  The  method  starts  with 
asking  the  experts  to  independently  arrive  at  a  forecast. 
Then  the  results  from  the  first  set  of  forecasts  are  sum¬ 
marized  as  a  simple  mean  and  range  (the  highest  and  lowest 
values)  (17:23).  Each  expert  is  then  asked  if,  based  on  the 
mean  and  range  of  the  previous  results,  he  or  she  desires  to 
keep  the  same  forecast  or  change  it.  Within  three  or  four 
iterations,  depending  on  the  nature  of  the  problem  and  the 
data,  a  consensus  is  arrived  at  and  the  experts  no  longer 
wish  to  change  their  forecasts  (17:23). 

In  this  study,  five  supply  officers  were  asked  to  act 
as  experts,  evaluate  the  available  data,  and  generate  values 
for  the  missing  data  points  using  their  knowledge  of  both 
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the  available  demand  data  and  forecasting  techniques.  At¬ 
tachment  4  lists  the  demand  data  used  in  the  exponential 
smoothing  forecasts,  while  Figures  10-21  offer  graphical 
representations  of  the  same  values.  Similarly.  Attachment  5 
lists  the  demand  data  used  in  the  Box-Jenkins  forecasts, 
while  Attachment  6  offers  a  graphical  representation  of  the 
same  values.  The  missing  data  points  replaced  by  Delphi¬ 
generated  values  are  underlined  in  Attachments  4  and  5. 
Recall  that  problems  with  obtaining  complete  data  forced  the 
use  of  slightly  different  data  sets  for  the  exponential 
smoothing  forecasts  and  the  Box-Jenkins  forecasts.  In 
addition,  recall  that  the  lack  of  complete  data  forced  the 
use  of  only  seven  National  Stock  Numbers  (NSN)  for  the  Box- 
Jenkins  forecasts. 

Research  Plan 

The  research  questions  were  addressed  in  the  following 
combination  of  consultation  with  experts,  literature  review, 
and  analysis  of  existing  data: 

1.  The  methods  by  which  the  SBSS  currently  addresses 
consumable  seasonal  demand  items  were  reported  in  the  liter¬ 
ature  review. 

2.  Other  methods  that  could  be  used  in  the  SBSS  to 
address  consumable  seasonal  demand  items  were  also  iden¬ 
tified  and  discussed  in  the  literature  review. 

3.  To  determine  if  items  currently  identified  as  sea¬ 
sonal  by  personnel  at  the  base-level  actually  display 
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seasonal  tendencies.  Stock  Control  personnel  at  Langley  AFB 
were  asked  for  a  list  of  items  that  are  given  special  man¬ 
agement  attention  based  on  suspected  seasonality.  The 
demand  histories  of  these  items  were  evaluated  to  determine 
if.  in  fact,  seasonal  trends  do  exist  in  the  data,  using  the 
autocorrelation  function. 

4.  To  determine  if  a  seasonally  adjusted  exponential 
smoothing  model,  a  simple  exponential  smoothing  model,  an 
AP.IMA  model,  or  the  current  SBSS  model  is  the  best  method 
for  forecasting  demand  for  consumable  seasonal  demand  items, 
the  demand  for  the  above  items  were  forecasted  using  these 
four  models  by  the  following  method. 

Demand  data,  obtained  from  AFLMC .  was  divided  into  two 
parts.  Three  years  of  data  was  used  to  define  the  para¬ 
meters  of  each  model  and  predict  the  consumption  of  the  item 
during  the  fourth  year,  which  was  used  to  test  the  models. 
The  best  model  will  have  the  lowest  Mean  Squared  Error 
defined  as  the  mean  of  the  squared  differences  between  each 
forecast  and  the  corresponding  actual  demand,  as  determined 
by  an  ANOVA  test  for  significant  differences. 

5.  In  order  to  identify  a  method  that  could  be  used  by 
base-level  personnel  to  identify  seasonal  items  one  maior 
difficulty  must  be  overcome.  Vital  to  any  system  of  iden¬ 
tifying  seasonal  data  are  demand  data  collected  over  three 
preferably  more  years.  Because  of  the  number  of  items 
handled  by  the  SBSS.  the  amount  of  demand  data  available  is 
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limited.  Therefore,  the  response  to  this  research  question 
is  necessarily  limited  to  recommendations  in  Chapter  V. 

Models 

SBSS  Forecasting  Methodology.  As  was  described  in 
Chapter  II.  the  daily  demand  rate  is  the  SBSS  demand  fore¬ 
cast.  Essentially  a  simple  moving  average  forecast,  it  is 
recalculated  as  needed  (26).  Whenever  a  demand  is  made  for 
an  item  in  the  SBSS.  an  R  is  placed  in  the  requirements 
computation  flag  field  of  the  item  record  (8:19-31).  If  a 
demand  level  has  already  been  calculated  that  quarter,  then 
requirements  computation  is  performed.  Requirements  com¬ 
putation  examines  how  many  of  the  item  are  available,  how 
many  are  needed,  and  whether  more  should  be  requisitioned  o 
some  should  be  returned  as  excess.  If  the  demand  level  has 
not  been  recalculated  for  the  item  that  quarter,  then  the 
demand  level  and  date  of  last  relevel  mg  would  be  updated 
(8:19-32).  Thus,  for  each  quarter,  the  demand  level  is 
calculated  only  once,  although  the  daily  demand  rate  is 
calculated  whenever  required  by  the  SBSS  (26). 

For  the  comparison  with  exponential  smoothing  and  Box- 
Jenkins.  the  DBR  was  calculated  as  follows.  The  SBSS  calcu 
lates  the  DDR  by  dividing  the  CRD  by  the  current  date  minus 
the  DOFD.  with  several  exceptions,  as  explained  in  Chapter 
(see  (11).  Two  adjustments  are  made  every  six  months.  The 
CRD  is  updated  to  the  current  DDR  times  the  lesser  of  the 
difference  between  the  current  date  and  DOFD.  or  365  (23:2) 
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So.  for  items  loaded  on  the  base  supply  system  for  more  than 
one  year,  the  old  CRD  is  changed  to  be  the  current  DDR  times 
365.  Also,  the  DOFD  is  adjusted  to  be  the  maximum  of  the 
DOFD  or  the  difference  between  the  current  date  and  365 
'23:2''.  Thus,  the  DOFD  is  adjusted  to  be  one  year  ago  for 
all  items  in  the  system  over  a  year. 

The  item  records  sent  by  AFLMC  were  from  the  end  of 
March  1989.  (Julian  date  9090).  Based  on  the  items'  DOFD  of 
8090  (365  days  prior  to  9090  day) .  the  item  records  had  been 
re  leveled  just  prior  to  being  transmitted  to  AFLMC.  The  CRD 
field  (obtained  from  the  item  record),  divided  by  365  equals 
the  DDR.  Also,  the  calculated  DDR.  multiplied  by  540.  gives 
the  CRD  after  the  reieveling  (new  CRD).  (According  to 
Patterson  540  days  are  the  maximum  number  of  days  in  a 
demand  history  (23:4).  although  this  is  not  exactly  equal  to 
365  +  180).  The  actual  demands  placed  on  each  item  were 
also  available  so  that  past  and  future  CRDs  could  be  ob¬ 
tained.  Thus,  the  DDRs  both  before  and  after  the  reieveling 
could  also  be  calculated  from  the  CRDs. 

A  sample  calculation  (for  NSN  1)  will  better  illustrate 
how  this  was  done.  Some  of  the  intermediate  and  final 
values  in  the  calculation  are  given  in  Table  III.  The  item 
record  (new)  CRD  was  divided  by  365  to  obtain  the  current 
DDR.  Since  the  DDR  before  the  reieveling  was  the  CRD  di¬ 
vided  by  540.  the  current  DDR.  .408.  was  multiplied  by  540 
to  obtain  220.  the  old  CRD.  Then  the  number  of  demands 
occurring  in  March  1989.  10.  was  subtracted  from  this  CRD  to 
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Table  III  DDR  Calculation  Values 


DEC  88 

JAN  89 

FEB  89 

MAR  89 

OLD  NEW 

APR  89 

CRD 

195 

203 

210 

220 

149 

154 

DDR 

0.433 

0.415 

0.413 

0.408 

0 . 408 

0 . 389 

ACTUAL 

8 

7 

10 

9 

9 

5 

DEMAND 

FCST 

13.43 

11.64 

12.8 

N/A 

12.24 

12.06 

obtain  210.  the  CRD  at  the  end  of  February  1989.  This  CRD 
was  divided  by  540  -  31  (the  number  of  days  of  demand  his¬ 
tory.  obtained  by  subtracting  the  number  of  days  in  March 
from  the  previous  number  of  days  demand  history)  to  obtain 
the  DDR  at  the  end  of  February.  Thus,  the  end  of  February 
DDR  was  multiplied  by  31  to  give  the  forecast  for  the  month 
of  March.  January  and  December  were  handled  the  same  way. 

The  calculation  of  the  SBSS  forecast  for  the  month  of 
April  was  more  straightforward.  For  the  April  forecast,  the 
new  CRD  was  divided  by  365  (the  new  current  date  -  DOFD)  and 
multiplied  by  30.  Although  not  used,  the  May  forecast  could 
be  obtained  by  adding  the  known  April  demands  to  the  new  CRD 
and  dividing  the  result  by  365  +  30  to  obtain  the  new  DDR. 
This  could  be  multiplied  by  31  to  obtain  the  May  forecast. 

Exponential  Smoothing  Methodology.  The  model  building 
for  the  Winters  and  simple  exponential  smoothing  models  was 
accomplished  using  the  commercially  available  MS-DOC  program 
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Forecast  Master.  The  program  uses  time  series  data  to 
create  models  which  it  uses  to  generate  forecasts. 

Of  the  four  exponential  smoothing  methods  available, 
two.  simple  exponential  smoothing  and  Winters  exponential 
smoothing  were  used.  Simple  exponential  smoothing  was  used 
as  being  somewhat  comparable  to  the  SBSS  method.  Both  meth¬ 
ods  produce  only  a  single  value  to  be  extended  for  all 
forecast  periods.  Winters'  seasonal  exponential  smoothing 
models  were  chosen  for  comparison  with  ARIMA  models  because 
both  compensate  for  seasonal  effects. 

The  technique  used  by  Forecast  Master  to  generate 
forecasts  is  quite  straightforward.  The  program  selects  the 
best  model  smoothing  parameters  as  determined  with  an  "iter¬ 
ative  search  method  to  minimize  squared  error  over  the 
historical  data"  (15:9-2).  Each  run  was  initialized  with 
the  same  smoothing  parameters  to  ensure  consistency.  For 
simple  exponential  smoothing  a  was  set  at  0.2.  For  Winters 
exponential  smoothing  a  was  also  set  at  0.2.  8  was  set  at 
0.05.  and  T  was  set  at  0.001.  After  the  computer  settled  on 
a  set  of  optimal  parameters  the  parameter  values  were  noted. 
The  forecast  generation  process  required  several  steps,  one 
for  each  forecast  .  The  first  forecast  was  for  the  first 
four  months  of  1989.  using  data  through  December  1988.  '  The 
second  forecast  was  for  February  through  April  1939  using 
data  through  January  1989.  and  so  on.  Thus,  four  separate 
sets  of  forecasts,  each  ending  in  April  1989  (the  end  of  th*3 
holdout  data),  were  made. 
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Box-Jenkins  Methodology.  The  three  major  steps  in¬ 
volved  in  the  Box-Jenkins  methodology  include  identifying  a 
tentative  model,  estimating  the  parameter  values  and  evalu¬ 
ating  the  model,  and  preparing  a  forecast. 

Identifying  a  Tentative  Model.  The  first  step  m 
building  a  Box-Jenkins  model  is  to  identify  a  tentative 
ARIMA  model.  As  Box  and  Jenkins  explain,  model  identifica¬ 
tion  is  an  inexact  science 

because  the  question  of  what  types  of  models  occur 
in  practice  and  in  what  circumstances,  is  a 
property  of  the  physical  world  and  cannot,  there¬ 
fore.  be  decided  bv  purely  mathematical  argument. 

(5:173) 

The  tentative  model  identified  is  only  to  be  used  if  diag¬ 
nostic  tests  demonstrate  the  model  fits  the  data  adequately. 

The  first  part  of  the  model  identification  process  is 
to  determine  if  differencing  is  necessary  to  make  the  data 
stationary.  Monstat ionar i ty  is  identified  by  autocorrela¬ 
tion  functions  that  do  not  die  off  (become  i ns igni f i cant ) 
but  appear  almost  linear  (5:174-175). 

The  second  part  of  the  process  is  to  determine  the  or¬ 
ders.  p  and  q.  of  the  AR  and  MA  processes .  Here  the 
autocorrelation  function  (ACF)  and  partial  autocorrelation 
functions  (PACF)  are  used.  Before  going  any  further,  how¬ 
ever,  the  partial  autocorrelation  function  must  be  des¬ 
cribed  . 

This  analysis  is  derived  from  that  given  by  Goodrich 
and  Stellwagon  in  their  guidebook  to  the  Scientific  Systems 
forecasting  program.  Forecast  Master  (15:8-9.10).  An  AR 
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process  may  be  identified  by  its  ACF .  although  p  remains 
unknown.  Since  an  AR  process  is  essentially  a  regression 
equation,  however,  p  may  be  determined  as  follows.  Assume  k 
to  be  the  order  of  the  regression  equation 


V  =  rr  v 

*  t  -H  1 

which  can  be 
In  addition. 


t-i  +  a2Yt-2  + 


written  and  the 
an  equation  with 


+  akYt-k 


values  of  the 
k+1  parameter 


+  et  [44] 

(15:8-10) 

as  determined. 

s  may  also  be 


written 


vt  aiYt-l 


+  a2 Yt-2  * 


+  a,  j.1 i  +  e^.  [45] 

"(15:8-10) 


If  a^  and  not  +  ^  is  significant  then  it  can  be  safely  as 
sumed  that  the  proper  order  of  the  equation  describing  the 
data  is  k.  The  PACF  is  a  function  describing  all  the  pos¬ 
sible  k  s.  (15:8-9. 10) 


If  the  autocorrelations  decline  exponentially  to  zero, 
then  an  AR  model  is  likely  to  best  explain  the  data.  Like¬ 
wise.  if  the  partial  autocorrelations  decline  exponentially 
to  zero,  then  an  MA  model  is  appropriate.  The  PACF.  then, 
will  be  zero  for  k  greater  than  p  and  non-zero  for  k  less 
than  p.  The  order  of  the  AR  model,  p.  is  determined  by  the 
number  of  "spikes"  or  significant  partial  autocorre lat ions . 
in  the  PACF.  For  example,  two  significant  spikes  indicate  a 
p  of  two.  (The  ACF  and  PACF  for  a  sample  AR  2  model,  m 
Figures  5  and  6.  show  significant  spikes  at  lags  1.  2.  3. 
and  6.  and  at  lags  1  and  2.  respectively.)  Similarly  the 
order  of  the  MA  model,  q.  is  determined  by  the  number  of 
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gure  6.  Partial  Autocorrelation  Function  for  Sample  AR  2  Series 


spikes  in  the  ACF:  two  significant  spikes  is  evidence  of  q 
equal  to  two.  (The  ACF  and  PACF  for  a  sample  MA  2  model,  in 
Figures  7  and  8.  show  significant  spikes  at  lags  1,  and  2. 
and  at  lags  1,  2.  3  and  4.  respectively.)  (5:175—177:  15:8- 
10:  16:62-65) 

Essentially  the  same  process  is  used  in  identifying 
seasonal  parameters,  except  the  ACF  and  PACF ' s  at  lags  12. 
24.  and  36  are  examined.  An  example  of  a  seasonal  pattern 
is  presented  in  Figure  9.  (Figure  9  shows  significant 
spikes  at  lags  1,  11.  12.  and  13.  with  the  spike  at  lag  12 
larger  than  those  at  lags  11  and  13).  The  number  of  data 
points  available  becomes  critical  in  this  area,  because  the 
identification  of  seasonality  depends  on  several  years  of 
data.  Also,  judgement  and  simply  trying  various  models 
become  important  analysis  tools.  (16:171-173) 

Parameter  Estimation  and  Model  Evaluation.  The 
next  step  in  the  Box-Jenkins  methodology  is  that  of  estima¬ 
ting  parameters  and  evaluating  the  model.  The  computer 
program.  Times,  performs  the  parameter  estimation  step  using 
a  nonlinear  (Marquardt)  routine  which  minimizes  both  "the 
relative  error  in  the  sum  of  squares"  and  the  relative  error 
in  each  parameter  value  (31:12). 

A  variety  of  diagnostic  tools  are  available  to  the  Box- 
Jenkins  forecaster.  The  first  of  these  is  to  check  for 
correlations  in  the  residuals.  If  the  differences  between 
the  actual  and  predicted  values  are  correlated  then  the 
model  is  not  capturing  all  the  underlying  patterns  (demand 
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gure  7.  Autocorrelation  Function  for  Sample  MA  2  Series 


PA c  Function  for  Sample  MA  2  Series 
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Figure  8.  Partial  Autocorrelation  Function  for  Sample  MA  2  Series 
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Figure  9.  Autocorrelation  Function  for  Sample  Seasonal  MA  and  Regular  MA  Series 


information).  Thus,  the  ACF  of  the  residuals  is  examined 
for  any  significant  spikes  (16:91-92). 

The  second  diagnostic  tool  is  to  plot  both  the  resid¬ 
uals  and  their  histogram.  The  residuals  should  have  a 
somewhat  norma  1  distribution,  be  random  (with  no  discernable 
pattern  about  the  mean)  and  have  a  mean  not  significantly 
different  from  zero.  According  to  Hoff,  the  last  can  be 
tested  by  determining  whether  or  not  the  residual  mean  is 
greater  than  "2  times  the  residual  standard  error  divided 
by  the  square  root  of  the  number  of  residuals"  (16:91) . 

A  third  diagnostic  tool  is  the  Q-statistic  (5:291-3). 
This  is  a  check  to  see  if  a  selected  group  of  res 'duals 
resulted  from  a  poor  model:  in  other  words,  whether  or  not 
they  are  correlated  as  a  group.  If  the  first  K  autocorrela¬ 
tions  are  used  in  the  formula 

K  2 

Q  =  n  2  ry{ a )  [ 46 ) 

k=l  (5:291) 

where  n  is  the  total  number  of  observations  minus  the  number 

2 

of  differences  taken  and  r^  is  the  residual  autocorrelation, 
and  Q  is  the  Q-statistic.  The  Q-statistic  is  compared  to  a 
chi-squared  distribution  with  (K-p-q)  degrees  of  freedom. 

If  the  Q  statistic  is  too  large  at  the  desired  confidence 
level,  then  the  model  may  be  regarded  as  inadequate.  (5:291- 
3) 

Another  diagnostic  tool  is  the  cumulative  periodogram 
(5:294-295).  This  is  a  means  of  evaluating  for  periodic 
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patterns  that  the  ACF  cannot  readily  detect.  The  periodo- 
gram  function  is  defined  as 

n  n 

Iff-)  =  2/n[(2  at,cos2irf  ■  t  )2  +  (2  atsin2irf  - 1 )  2  [47] 

t=l  '  *  t=l  (5 : 294) 

where  fi  =  i/n  is  the  frequency.  If  the  series  in  question 

is  random,  then  the  periodogram  takes  on  the  appearance  of  a 

straight  line  running  from  the  coordinates  (0.0)  to  (0.5,0). 

A  sufficiently  large  deviation  from  this  line  is  indicative 

of  periodicity  not  captured  by  the  model. 

A  model  passing  these  basic  tests  may  be  considered 
adequate  for  the  forecaster  to  move  on  to  the  final  step: 
using  the  model  to  prepare  forecasts.  The  process  of 
generating  forecasts  is  relatively  straightforward.  Once 
the  parameters  and  form  of  the  forecast  model  are  known,  the 
computer  simply  solves  the  equation  to  obtain  the  forecast. 

Forecast  Comparisons.  Once  forecasts  were  obtained 
from  each  forecast  method,  they  were  compared  with  the 
holdout  sample  from  the  f-irst  four  months  of  1989  (January 
1989-April  1989).  Each  model  was  evaluated  by  MSE  on  how 
well  it  predicted  the  values  in  the  holdout  sample.  The 
forecasts  from  each  model  were  obtained  four  times  and  the 
number  of  times  the  model  had  the  lowest  MSE  was  recorded  to 
indicate  how  successful  each  technique  was  at  predicting  the 
values  in  the  holdout  sample.  The  first  forecast  used 
demand  data  through  December  1988  to  forecast  demands  for 
January  through  April  1989.  The  second  forecast  used  demand 
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data  through  January  1989  to  forecast  demands  for  February 
through  April  1989  and  so  on.  Since  only  seven  NSNs  could 
be  used  to  create  ARIMA  models  because  of  the  two  year  data 
gap  problem,  only  those  seven  NSNs  were  used  to  compare  all 
four  models.  The  other  five  NSNs  were  not  used  in  this 
evaluation  step. 

Summary 

This  chapter  discussed  the  data  acquisition  and  manip¬ 
ulation  process.  It  also  highlighted  the  methods  used  to 
actually  make  forecasts  from  each  forecasting  method  and  how 
those  forecasts  would  be  compared  against  the  holdout  sam¬ 
ple.  The  next  chapter  continues  with  the  results  of  the 
forecasts  and  how  each  model  fared  in  predicting  the  holdout 
samp  1 e . 
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Introduct ion 


This  chapter  discusses  the  results  of  the  data  analysis 
and  the  relationships  found  during  the  study.  First  results 
of  the  seasonal  tests  will  be  discussed.  Next,  the  exponen¬ 
tial  smoothing  models  will  be  discussed,  followed  by  the 
ARIMA  models.  The  relationships  between  the  models  of 
similar  stock  numbers  will  also  be  discussed.  This  will  be 
followed  by  the  results  of  the  forecasts  and  the  rating  of 
each  model  on  its  performance  in  predicting  the  holdout 
data . 

Seasonality  Tests 

Two  methods  of  determining  whether  data  is  seasonal  or 
not  were  mentioned  in  the  previous  chapter:  examination  of 
plots  and  evaluation  of  the  autocorrelation  function. 

Demand  Data  Evaluation.  First,  data  plots  will  be 
evaluated.  Recall  from  Chapter  III  that  the  data  used  to 
build  the  exponential  smoothing  models  is  slightly  different 
from  that  used  to  build  the  ARIMA  models,  because  of  the 
problems  in  obtaining  enough  data.  Since  ARIMA  models  were 
built  for  only  seven  of  the  twelve  items,  this  section  on 
analysis  of  seasonality  examined  the  demand  data  used  in 
building  exponential  smoothing  models,  so  that  all  twelve 
items'  demands  could  be  evaluated.  The  exponential  smooth¬ 
ing  demand  data  are  plotted  in  Figure  10  through  Figure  21. 
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demand  data 


DEMAND  DATA  FOR  EXPONENTIAL  SMOOTHING 
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Exponential  Smoothing  Forecast  Demand  Data  for  NS 
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Exponential  Smoothing  Forecast  Demand  Data  for  NSN 
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l5.  Exponential  Smoothing 
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Jinoo thing  Forecast  Demand 


DEMAND  DATA  FOR  EXPONENTIAL  SMOOTHING 


SQNW0O 
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Figure  17.  Exponential  Smoothing  Forecast  Demand  Data  for  NSN 


1  NG 


demand  DATA 


DEMAND  DATA  FOR  EXPONENTIAL  SMOO 


SONW0O 


re  20.  Exponential  Smoothing  Forecast  Demand  Data  for  NSN  11 


Attachment  5  lists  the  data  used  to  create  the  plots.  For 
the  sake  of  completeness,  the  demand  data  used  in  the  Box- 
Jenkins  methodology  to  generate  the  ARIMA  models  are  listed 
in  Attachment  6  with  the  associated  graphs  in  Attachment  7. 

Examination  of  the  demands  used  to  build  the  exponen¬ 
tial  smoothing  models  reveals  some  consistency  in  which 
months  had  the  highest  demands.  Table  IV  summarises  which 
months  had  the  high  demands  for  each  item.  The  months  or 
September.  October  and  February  had  the  largest  number  of 
high  demands.  February  1987  had  the  most  high  item  demands: 
MSN ' s  4.  11.  and  12  had  their  highest  demands  in  that  month. 
Three  months  had  the  highest  demands  for  two  items:  Septem¬ 
ber  1986  had  more  demands  than  the  other  months  for  NSNs  1 
and  9;  February  1^88  had  more  demands  than  the  other  months 
for  NSNs  5  and  7 :  October  1988  had  more  demands  than  the 
other  months  for  NSN ' s  2  and  10.  Three  other  months  had  the 
highest  demands  for  the  final  three  items:  demand  for  NSN  6 
peaked  in  October  1987.  demand  for  NSN  3  peaked  in  November 
1987.  and  demand  for  NSN  0  peaked  in  January  1989.  (Dis¬ 
counting  that  value  as  a  holdout  value  would  move  the  peak 
demand  for  NSN  8  to  November  1907.)  In  short,  the  fall  and 
winter  months  experienced  a  much  greater  demand  for  these 
i terns . 

The  following  discussion  examines  the  demand  patterns 
of  each  NSN.  As  was  noted  in  Chapter  ITT.  plotted  seasonal 
data  has  simi lar  shaped  demand  curves  from  year  to  year  and 
is  consistent  in  increasing  and  decreasing  from  year  to 
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Table  IV.  NSN  Versus  Period  of  Highest  Demands 


NSN  NOMENCLATURE  PERIOD  OF  PEAK  DEMANDS 

HIGHEST  2ND  HIGHEST 


NSN'i 

BOOT.  MUKLUK  LARGE 

SEP  86 

FEB 

88 

NSN  2 

PARKA. ECW  GRN  X-SMA 

OCT  88 

FEB 

87 

NSN  3 

SOCKS. MENS  C/W  NA10 

MOV  87 

OCT 

87 

NSM4 

COAT  MANS  SML-SHOR 

FEB  87 

JAN 

87 

NSN5 

COAT  MAN  MED  SHT 

FEB  87 

NOV 

87 

NSN6 

SOCKS. MENS  C/W  NA14 

OCT  87 

SEP 

88 

NSN7 

COAT  MAN  MED  LG  33 

FEB  88 

JAN 

88  : 

NSN8 

UNDERSHIRT  EX  WE A  L 

JAN  89 

NOV 

87 

NSN9 

GLOVE  SHELLS.  CW  S22 

SEP  86 

OCT 

88 

NSN  10 

SOCKS . MEN . CW . NAUTR9 

OCT  88 

NOV 

87 

NSN  11 

BOOT  MUKLUK  MANS  SM 

FEB  87 

OCT 

87  < 

NSN  12 

BOOT  MUKLUK  MED 

FEB  87 

FEB 

88 

year.  An  overlay  of  a  plot  of  data  from  one  year  to  the 
next  would  show  little  or  consistent  change. 

NSN  1  had  its  peak  demand  of  55  in  September  1986. 
followed  by  the  second  highest  demand  of  47  in  February 
1988.  The  next  highest  demands  of  39  and  37  were  in  January 
and  February  1987.  respectively,  followed  by  the  5th  highest 
demand  of  32  in  October  1987.  This  pattern  of  dual  peak 
demands  continued  in  1988  and  1989  with  the  highest  demands 
(only  27  and  10.  respectively)  again  in  September  and  Febru¬ 
ary.  The  lowest  demands  of  0  were  in  June  and  July  of  both 
1987  and  1988.  and  March  and  April  of  1988.  However,  in 
1987  over  30  items  were  demanded  in  March  and  April.  Thus 
NSN  1  had  a  fairly  well-established  pattern  of  dual  peak 
demands  with  few  demands  in  June  and  July. 

NSN  2  has  a  similar'  demand  pattern,  but  with  a  smaller 
number  of  demands,  wit-h  the  peak  demand  of  9  in  October 
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1988.  followed  by  the  2nd  highest  monthly  demand  of  8  in 
February  1987.  During  the  three  months  of  October  1987.  and 
March  and  April  1988.  the  monthly  demands  were  6.  As  with 
NSN  1.  however  there  were  relative  peak  demands  of  3  in 
November  1987  and  March  iu86.  followed  by  an  October  1988 
peak.  With  such  low  peak  demands,  no  real  pattern  of  low 
demands  could  be  found.  So.  the  pattern  of  dual  peaks, 
several  months  apart,  is  also  evident  in  NSN  2's  demands. 

NSN  3  has  a  more  variable  demand  pattern.  The  highest, 
demand  of  533  occurred  in  November  1987.  with  the  second 
highest  demand  of  325  occurring  in  October  1987.  The  third 
and  fifth  highest  demands  of  236  and  154.  respectively,  also 
occurred  in  consecutive  months:  September  and  October  of 
1988.  To  continue  the  pattern,  the  fourth  and  sixth  highest 
demands  were  also  in  consecutive  months,  with  21 J  demanded 
in  February  1987.  and  153  demanded  in  January.  The  varia¬ 
tion  in  pattern  primarily  occurred  in  late  1986.  with  rel¬ 
atively  high  demands  ror  four  consecutive  months.  This 
plateau  in  demand  is  unmatched  in  other  years.  Lev;  demands, 
those  in  the  single  digits,  were  consistently  in  the  months 
of  June  and  July.  In  all  three  years,  then,  demands  were 
high  in  October  and  November,  but  in  1986  and  early  198"7. 
the  peak  demands  for  the  year  were  spread  out  and  the  high¬ 
est  peak  demands  delayed. 

For  NSN  4.  the  demands  followed  d^cr^asing  pattern  ov<ir 
the  years.  The  winter  (defined  as  th**  period  of  r^ak  de¬ 
mands  occurring  roughly  from  Sept  *»mb»r  ‘-f  y«ar  through 
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March  of  the  following  yean  of  1986  and  1987  had  over  55% 
of  the  demands  for  the  period.  All  of  the  highest  demands 
for  the  item  occurred  between  September  1986  and  April  1987. 
with  the  highest  demand  of  44  in  February  1987.  The  next 
’■wo  hign^st  demands  were  for  86  items  in  January  19*?“.  and 
33  items  in  December  1986.  September  1986  had  the  next 
highest  demand  at  30.  while  November  1986.  and  March  and 
April  1987  had  22  demands  each.  This  block  of  higher  de¬ 
mands  was  somewhat  repeated  a  year  later  with  the  high 
demand  for  the  year  of  18  again  occurring  in  February  1986. 
For  the  winter  of  1968-1989.  however,  the  pattern  was  not 
repeated.  The  peak  demands  of  13  occurred  in  September  and 
October,  and  demands  declined  fairly  steadily  after  that 
point.  Again,  low  demands  were  very  consistent,  occurring 
in  the  months  of  May.  June,  and  July,  although  November 
1988.  and  March  and  April  1989  also  had  demands  of  only  1. 
Thus,  for  MSN  4.  declining  demands  followed  a  pattern  of 
plateaumg  for  the  winter  months  for  two  of  the  three  years. 

MSN  5 '  s  demands,  follow  somewhat  a  combination  of  the 
pattern  of  NS Ns  1  and  2.  The  peak  demand  of  139  occurred  in 
February  .  with  the  second  highest  peak  demand  of  107  m 

November  1987.  Much  lower  demands  of  58  in  February  L98~ 
and  38  in  January  1987  are  the  third  and  fourth  highest 
demands  of  the  previous  winter  period.  The  fifth  highest 
demand  of  31  occurred  in  both  September  and  October  T^t 
Lew  demands  for  N3N  5  were  somewhat  consistent.  They  :■<: 
rurred  from  May  through  September  in  1987  and  March  thrown!; 
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August  in  1988,  although  other  months,  such  as  February  and 
April  1989  also  had  low  demands.  So.  again,  there  is  a 
pattern  of  high  demands  in  the  fall,  a  drop  in  demands 
throu  ;h  December  and  January,  followed  by  a  high  demand  in 
February  as  in  NSNs  1  and  2's  demands.  This  pattern  ends  in 
the  t’nal  winter  of  demand  data  with  two  relative  peak 
demands  of  10  and  13  in  September  and  November  1988  respec¬ 
tively.  but  with  no  follow-on  peak  demand  in  February.  As- 
in  th  case  of  NSNs  1  and  4.  the  demands  almost  established 
a  pattern  for  the  first  two  years  which  changed  for  the 
wi :  ter  of  1988-1989. 

MSN  6  followed  even  less  of  a  pattern.  The  peak  demand 
of  93  was  in  October  1987.  while  the  second  highest  demand 
of  57  was  in  September  '988.  The  third  highest  demand  of  29 
occurred  in  December  1988.  while  the  fourth  highest  demand 
of  21  was  in  February  1:87.  The  fifth  highest  demand  was  20 
in  November  1987.  In  the  winter  of  1986-1987  the  peak 
demand  was'  in  February.  The  following  winter  (1987-19881 
the  peak  demands  were  in  consecutive  months.  October  and 
November.  The  final  winter  the  peak  demands  were  split  into 
th“  two  months  of  September  and  December.  The  low  demands 
are  also  inconsistent.  In  1987.  they  extend  from  April 
through  September,  whil  in  1988  they  extend  from  December 
f',987)  through  July.  3<  •  there  is  a  less  clear  pattern  to 
the  demands  of  MSN  6. 

MSN  7  al  :o  failed  t-  show  a  clear  demand  pattern.  The 
peak  de  and  f  90  occurred  in  February  L9S8.  The  second 


81 


highest  demands  of  30  occurred  during  January  1988  and 

Sepi' ember  L986 .  The  fourth  highest  demands  of  18  occurred 

In  : •ecember  1966.  and  January  and  February  1987,  The  next 

hig.nest  d  emand  kc  event  hi  occurred  in  Serf  ember  1988.  A 

mlx-'d.  "•mix  paf'.ern  might  be  cr leaned  from  these  demand  data  » 

if  the  highest  peak  of  over  90  were  not  so  high.  Peak 

demands  recurred  twice  in  September,  and  twice  in  January 

and  February.  The  low  demands  are  relatively  consistent . 

running  from  March  through  September  in  1967  and  March 

through  August  in  1968.  In  conclusion,  there  is  not  much  o-f 

a  visual  pattern  to  the  NSN  7  demands. 

MSN  a  had  by  far  t he  highest  demand  rate  of  any  of  the 
MS!!'-  s *  : .i d  i  e d .  but  little  p a tter n  .  .The  tri o n thly  d e m a n d s 
ranged  -om  2  to  873.  the  highest  demand  occurring  in  Janu¬ 
ary  iG8-  The  second  highest  demand  of  390  fell  in  November 
1987,  while  the  third  highest  demand  of  502  occurred  :n 
December  1989.  October  1988  had  the  fourth  highest  demand 
of  412  and  February  1987  had  the  fifth  highest  demand  of 
350.  The  lowest  monthly  demands  were  fairly  consistent,  o -- 
earring  in  May  and  June  of  1987  and  June  and  July  of  I960. 

Other  than  this,  no  obvious  pattern  reveals  itself  in  the 
da'  -i.  The  winter  1366-1987  demand  data  consists  of  a  number 
of  relatively  high  demands  from  September  to  March,  the  two  ' 

highest  occurring  in  September  and  February .  The  following 

* 

•winter's  d e m a n d  data  is  c h a r a c t c r  i  r. e d  by  a  1  a i  g e  p e a k  in 
Mov'-mber  amidst  a  saw-tooth  pattern  with  relatively  large 
demands  in  Angus'-.  October.  January,  and  March.  The  winter 
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of  1988  and  1989  Is  still  different  with  two  large  peaks, 
one  centered  around  October,  the  other  centered  around 
January.  So,  in  no  two  years  are  the  peak  demands  in  the 
same  months,  nor  are  the  patterns  of  demand  similar  in  any 
'•f  the  three  winters 

For  N3N  9.  the  demand  pattern  is  more  apparent.  The 
peak  demand  of  92  occurred  in  September  1986.  The  second 
highest  demand  of  72  occurred  in  October  1988.  while  the 
third  highest  demand  of  56  occurred  in  September  of  the  same 
year.  The  fourth  highest  demand  of  52  followed  in  January 
1989.  while  the  fifth  highest  demand  of  27  took  place  in 
January  1987.  The  lowest  demands  were  not  consistent.  In 
1987  the  lowest  demands  occurred  only  in  June.  July,  and 
August  (and  November),  while  in  1988  February  through  August 
had  low.  demands.  Smaller,  but  relatively  high  demands  also 
occurred  in  December  of  i‘^87  and  January  of  1988  along  with 
smaller  yet  relatively  high  demands  in  September  and  October 
of  1987.  The  pattern  is  fairly  evident  although  not  consis¬ 
tent  in  magnitude,  since  the  winter  1987-1988  peaks  are  so 
low.  For  NSN  9.  the  high  demands  occur  in  January  fairly 
consistently,  with  another  peak  in  demands  in  September  or 
October  or  both  September  and  October. 

The  NSN  10  demand  pattern  is  similar  to  that  of  NSN  9. 
The  highest  demand  of  277  occurred  during  October  1988. 
while  the  second  highest  demand  of  206  occurred  in  November 
1987.  The  third,  fourth  and  fifth  highest  demands  of  129. 
107  and  105  respectively  occurred  in  February  1987. 
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September  1986.  and  January  1987.  The  lowest  demands  oc¬ 
curred  in  May  through  September  1987  (with  the  exception  of 
a  higher  August  value),  while  occurring  from  only  April 
through  July  of  1988.  The  high  demands  were  mere  consis¬ 
tent.  occurring  m  October  and  November,  and  to  a  lesser 
extent  in  January  or  February.  Each  winter  period,  espe¬ 
cially  1986-1987  and  1988-1989.  had  at  least  a  double  peak 
demand.  The  winter  1987-1988  demands  showed  a  pattern  of 
one  central  peak  demand  in  November .  a  lesser  one  in  August 
and  another  lesser  one  in  March. 

The  demand  data  for  NSN  11  has  a  less  clear-cut  pat¬ 
tern.  The  highest,  demand  of  37  took  place  in  February  1987. 
The  next  highest  demand  of  30  occurred  in  both  October  198"7 
and  October  1988.  The  fourth  highest  demand  of  28  occurred 
in  May  1988.  while  the  fifth  highest  demand  of  27  occurred 
in  September  1988.  The  low  demands  for  198'7  took  place 
April  through  September,  but  in  the  following  year  occurred 
during  February  through  August,  with  the  previously  men¬ 
tioned  exception  of  May  1988.  The  winter  1986-1987  demand 
data  have  dual  peaks,  one  in  September,  the  other  in  Febru¬ 
ary.  The  next  winter,  dual  peak  demands  occurred,  one  m 
October,  the  second  in  January.  These  were  followed  by  an 
unusual  peak  in  May  of  1988.  which  was  subsequently  followed 
by  still  two  more  peaks  in  September  and  October  and  a  minor 
one  in  January  1989.  Excepting  the  May  peak.  September-. 
October,  and  January  tended  to  have  peak  demands,  but  there 
is  no  apparent  consistency  to  their  relative  heights  because 


of  the  changes  in  which  months  had  the  greatest  demands  from 
year  to  year . 

NSN  12  has  a  similar  demand  pattern.  The  highest 
monthly  demand  of  79  took  place  in  February  1987.  The  next 
highest  demand  of  68  occurred  in  February  1988.  followed  by 
the  third  highest  demand  of  66  in  September  1986.  The 
fourth  highest  demand  of  65  occurred  in  October  1988  and  the 
fifth  highest  demand  of  68  occurred  in  January  1987.  The 
low  demands  occurred  April  through  September  in  1987  (with 
August  an  exception)  and  -in  April  through  August  in  1988. 
Again  the  pattern  of  dual  peaks  occurring  inconsistently 
throughout  the  year  is  evident.  One  peak  occurs  m  Septem¬ 
ber  or  October  followed  by  another  in  January  or  February, 
but  there  is  no  consistency  tp  the  relative  heights  of  the 
peaks  for  given  months  from  year  to  year. 

G e ne r a  1  Dema n d  Patter n s  .  The  most  striking  aspect 
of  the  plots  ar^  the  peaks  which  generally  occur  in  the  fall 
and  winter  months,  followed  by  low  demands  in  the  spring  and 
summer  months.  A  pattern,  although  sometimes  slight,  exists 
in  the  demands  of  NSNs  1.  2.  3.  5.  9.  10.  11.  and  12.  (Note 
that  these  are  the  same  NSNs.  with  the  addition  of  NSN  2. 
that  showed  enough  of  a  consistent  pattern  to  "pass"  the 
Wilcoxon  Signed  Rank  Test.)  The  pattern  consists  of  high 
demand  characterised  by  a  peak  during  one  or  more  .-f  the 
months  September  through  November,  followed  by  low  demands 
in  December  and  another  peak  in  January.  February,  or  March. 
NSN  4 ’ s  demands  had  a  pattern  of  plateaus  instead  of  peaks 
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during  the  months  December  through  February.  The  other 
NSNs .  6.  7.  and  8  had  little  or  no  obvious  pattern. 

Several  possible  reasons  exist  for  the  demand  data 
differing  from  year  to  year.  First,  organizations  tend  to 
buy  odd  weather  gear  in  large  quantities  only  a  few  t  linos  i 
year,  so  year  to  year  changes  in  when  large  organizations 
bought  winter  gear  could  easily  change  monthly  demands  from 
year  to  year.  Second,  the  supply  personnel  could  change  the 
month  of  peak  demands  from  year  to  the  next  by  changing  when 
they  ask  organizations  to  request  winter  gear,  and  by  order 
ing  differing  quantities  according  to  what  the  response  is. 

A  third  possibility  is  that  changing  budgets  could  cause 
changes  in  how  the  organizations  choose  the  spend  their 
money.  Years  with  low  demands  across  the  spectrum  of  -old 
weather  gear  could  mean  restrictions  on  who  was  authoriz'd 
to  possess  an  item  and  how  many  of  a  given  item  an  organiza¬ 
tion  could  procure,  or  a  large  decrease  m  the  manning  <>f 
the  organization.  A  different  possibility  is  that  nuid 
weather  one  year  could  affect  demands  both  that  year  and  the 
next.  Fewer  items  will  be  used  during  a  mild  winter  and 
those  that  are  used  will  be  used  less  frequently  so  replace¬ 
ments  would  be  fewer  the  following  year.  Another  possible 
cause  for  changing  demand  is  replacement  of  the  item  by 
another,  either  because  it  falls  out  of  fav-*r.  or  because  i* 
is  simply  replaced  by  another  NSN .  This,  however,  is  un¬ 
likely  for  the  items  considered  here .  ftill  another  pos¬ 
sibility  is  that  the  extremely  low  demand  rates  of  some 
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items  just  do  not  provide  sufficient  data  to  make  an  under¬ 
lying  demand  pattern  visible. 

Grouped  Item  Demand  Evaluation.  In  addition  to 
the  individual  items,  three  different  sets  of  grouped,  like, 
items  were  evaluated.  Thr^e  sices  each  of  mukluks.  cc  I'd 
weather  seeks,  and  coats  were  among  the  items  selected  for 
ana  lysis . 

The  mukluk  demands,  cold  weather  sock  demands,  and  coat 
demands  are  presented  in  Figure  23.  Figure  24.  and 
Figure  22.  respectively.  An  examination  of  the  demand  pat¬ 
terns  for  like  items  showed  some  similarities  and  some 
differences.  The  mukluks  (NSNs  1.  11.  and  12)  were  repre  - 
sented  by  the  same  basic  demand  pattern  of  two  or  three 
peaks  of  varying  height  over  each  winter’s  heavy  demand 
period . 

The  socks'  (NSNs  3.  b.  and  10)  demand  patterns  were 
more  similar  over  the  three  winter  period.  The  first  winter 
had  a  plateau  shaped  demand  pattern  followed  by  a  single 
large  peak  the  next  winter,  and  finally  a  single  peak  the 
third  winter  followed  by  relatively  level  demands  or  a 
single  additional  peak.  The  coats’  (NSNs  4.  5.  and  ?) 
demand  patterns  were  less  similar.  NSNs  5  and  7  had  basi¬ 
cally  the  same  demand  pattern  with  plateaus  and  peak 
switched  during  the  first  winter,  and  two  peaks  in  ♦-h*a 
second  winter’s  demand  of  NSN  5  versus  only  one  for  NSN 
However.  NSN  4  followed  an  slightly  different  demand  pat  ♦■air. 
with  no  peak  demand  during  the  second  winter.  So.  for  each 
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of  the  three  groups,  although  some  similarities  could  be 
found  in  the  demand  patterns  of  like  items,  the  similarities 
were  limited. 

Autocorrelation  Function  Evaluation.  The  limited 
similarities  in  the  demand  patterns  of  like  items  are  also 
evident  in  the  ACFs  of  the  items'  demand  data  shown  in 
Attachment  7.  As  an  example,  the  ACF  of  NSN  1  is  displayed 
in  Figure  25.  The  calculated  ACFs  are  presented  in  a  bar 
rormat ,  while  the  upper  and  lower  bounds  of  significance  for 
the  ACFs  at  each  lag  are  presented  as  lines.  These  sig¬ 
nificance  bounds  are  set  according  to  the  rule  that  for  the 
first  three  lags  an  autocorrelation  greater  than  one-half 
the  standard  error  is  significant,  for  the  next  three  lags 
an  autocorrelation  greater  than  the  standard  error  is  sig¬ 
nificant.  for  lags  7-12  an  autocorrelat ion  greater  than 
twice  the  standard  error  is  significant,  and  for  higher 
numbered  lags  an  autocorrelation  greater  than  three  times 
the  standard  error  is  significant  (7). 

None  of  the  NSNs  had  a  strong  enough  seasonal  demand 
pattern  to  generate  a  statistically  significant  seasonal 
ACF.  although  some  were  close  to  being  significant.  For 
example,  NSNs  5  and  8  had  relatively  strong  autocorrela¬ 
tions  at  lag  12,  but  in  both  cases  they  were  not  statisti¬ 
cally  significant.  Interestingly,  NSNs  10  and  12  had  high 
autocorrelations  at  the  11th  lag.  This  implies  a  demand 
pattern  shifting  from  one  month  to  the  next  each  year. 
Finally,  NSN  2  had  a  demand  pattern  with  a  relatively 
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Autocorrelation  Function  for  NSN 


significant  5th  and  6th  lag.  leading  to  the  conclusion  that 
demands  were  on  a  six  month  cycle  instead  of  the  expected  12 
month  cycle. 

The  net  results  of  the  seasonal  tests  were  that  the 
demand  data  are  not  strongly  seasonal.  Although  the  demand 
data  appeared  superficially  seasonal,  closer  analysis  showed 
the  demand  data  generally  do  have  some  seasonal  patterns, 
although  they  are  not  very  strong.  The  evaluation  of  the 
ACFs  gave  similar  results.  None  of  the  demand  data  had  a 
strong  seasonal  pattern.  Having  examined  the  seasonality  of 
the  items,  the  performance  of  the  forecast  models  will  be 
discussed  next. 

Exponential  Smoothing  Models 

Two  types  of  exponential  smoothing  models  had  their 
forecasting  abilities  tested.  The  first  was  simple 
exponential  smoothing  [14] .  The  second  was  Winters'  sea¬ 
sonal  exponential  smoothing  [211.  [22].  [231.  and  f 24 1 .  For 

each  item  a  value  for  the  parameters  (a  for  simple  exponen¬ 
tial  smoothing:  and  a.  (3,  and  T  for  seasonal  smoothing)  w^r-' 
generated.  These  values  are  listed  in  Table  V. 

A  further  analysis  was  performed  on  the  three  groups  of 
like  items  to  determine  if  any  relationships  could  be  found 
between  the  exponential  smoothing  models  for  each  item  in 
each  group.  To  determine  if  the  parameters  for  the  exponen¬ 
tial  smoothing  models  were  the  same  for  th-  different  sixes 
within  each  group  a  one  sample  t-test  was  performed  on  the 
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Table  V .  Exponential  Smoothing  Parameters 


EXPONENTIAL  SMOOTHING  PARAMETERS 
SINGLE  WINTERS' 


NSN 

a 

a 

B 

r 

1 

X 

.  028 

.  050 

.  046 

.  332 

n 

.  044 

.  228 

.  028 

.410 

j 

.018 

.  030 

.  030 

.  244 

4 

.  684 

.  073 

.  061 

.  391 

5 

.  038 

.  038 

.  000 

.  251 

6 

.  054 

.  040 

.  055 

.  075 

7 

.  035 

.  035 

.  026 

.  2^7 

8 

.  076 

.  061 

.  050 

.  229 

9 

.004 

.  000 

.  033 

.  266 

10 

.044 

.  027 

.  003 

.163 

11 

.  031 

.012 

.013 

.  271 

12 

.  848 

.  062 

.  020 

.  367 

simple  exponential  smoothing  parameter,  while  a  one  way 
analysis  of  variance  (AOV)  was  performed  on  the  three 
parameter  values  for  the  seasonal  model.  The  results  of 
those  tests  are  in  Table  VI. 

For  the  mukluks  the  1.14  t-test  value  for  the  simple 
exponential  smoothing  models  indicates  that  the  three  simple 


exponential  smoothing  models  are  not  statistically  dif¬ 
ferent.  The  large  F  value  of  75. 55.  however,  leads  to  the 
conclusion  that  the  Winters'  seasonal  exponential  smoothing 
models  for  the  three  mukluk  sixes  are  indeed  different. 

For  the  socks.  (NSNs  3.  6.  and  10).  the  results  ai  e 
different.  The  simple  exponential  smoothing  models  are 
st a t ist l ca 1 ly  different  at  the  05 %  confidence  level.  The 
Winters'  seasonal  exponential  smoothing  models,  however,  are 
statist » ca  1  l  y  different,  at  a  95%  confidence  level. 
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Table  VI. 
Models 

T-test 

and  AOV  Results  for 

Exponential  Smoothing 

GROUP 

NS  Ns 

EXPONENTIAL 

SIMPLE 

SMOOTHING 

WINTERS ' 

T-VALUE 

F-VALUE 

MUKLUKS 

1.11.12 

1  .  14 

75.55 

SOCKS 

3.6.10 

3.60 

6.43 

COATS 

4,5,7 

1  .  17 

31.11 

In  the  case  of  the  coats,  their  simple  exponential 
models  are  not  statistically  different  at  a  95%  confidence 
level.  However,  one  could  be  confident  in  saying  the  three 
items  could  not  be  represented  by  a  single  Winters'  season¬ 
al  exponential  smoothing  model  at  the  95%  confidence  level. 

Given  the  above  information,  then,  a  base  could  choose 
to  use  a  single  simple  exponential  smoothing  model  to  fore¬ 
cast  demands  for  mukluks.  another  model  for  coats,  but  would 
need  three  models  for  socks.  Since  the  sample  is  so  small, 
one  cannot  tell  if  this  would  hold  true  for  all  groups  of 
similar  but  different  sized  items.  It  must  be  noted  that 


the  accuracy  of  these  models  has  not  yet  been  evaluated. 

In  the  case  of  Winters'  seasonal  exponential  smoothing, 
one  model  could  not  be  considered  appropriate  for  the  dif¬ 
ferent  sizes  of  mukluks,  socks,  and  coats.  It  is  reasonable 
that  one  Winters'  seasonal  exponential  smoothing  model  is 
less  likely  to  work  for  more  than  one  item  since  there  are 
three  parameters  to  match,  not  just  one. 
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The  net  result  of  the  exponential  model  analysis  is 
that,  based  on  the  limited  sample  used,  one  simple  exponen¬ 
tial  smoothing  model  may  be  used,  in  some  cases,  to  forecast 
different  sizes  of  similar  items,  but  not  one  Winters' 
seasonal  exponential  smoothing  model.  For  each  NSN  a  sep¬ 
arate  Winters'  seasonal  exponential  smoothing  model  was 
required.  One  possible  reason  is  that  with  the  difference 
in  demand  patterns  between  similar  items  identified  earlier, 
the  more  sensitive  Winters'  models  followed  demands  for  each 
item  so  closely  they  could  not  reasonably  model  other  items' 
d emands . 

Next.  ARIMA  models  will  be  considered  to  observe  how 
they  model  the  demands  for  seven  of  the  sample  items.  Each 
ARIMA  model  will  be  compared  with  the  forecasts  generated  by 
the  simple  and  Winters'  seasonal  exponential  smoothing 
mode  1 s . 

ARIMA  Models 

Only  seven  items  were  forecast  using  ARIMA  models 
because  of  the  break  in  the  data  for  the  period  1984  through 
the  end  of  1986.  Table  VII  contains  a  summary  of  the  models 
obtained  through  the  Box-Jenkins  methodology,  while  Attach¬ 
ment  8  lists  the  diagnostic  outputs  from  the  computer. 
Recall,  AR  and  MA  stand  for  autoregressive  and  moving  aver¬ 
age  respectively,  while  RDIF  stands  for  regular  differenced. 
No  twelve  period  seasonal  parameters  were  necessary  to 
the  data,  nor  was  any  seasonal  differencing  necessary .  with 
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Table  VIT.  Box— Jenkins  Methodology  Models 


PARAMETER 

TYPE 

AND  VALUE 

NSN 

1 

AR  1 

-0.53133 

(A) 

AR  2 

-0 . 36419 

RDIF 

1 

NSN 

2 

N/A 

NSN 

3 

MA  1 

-0.28383 

( B) 

MEAN 

104.38 

NSN 

4 

N/A 

(C) 

NSN 

5 

AR  3 

0.62381 

(C) 

AR  6 

-0.47797 

MA  1 

0.42724 

MA  2 

0 . 26339 

RDIF 

1 

NSN 

6 

N/A 

(B) 

NSN 

7 

N/A 

(C) 

NSN 

8 

N/A 

MSN 

9 

AR  1 

0 . 38984 

AR  3 

0.26171 

MA  2 

0 . 24589 

MEAN 

20.036 

NSN 

10 

MA  1 

0.65234 

(B) 

MA  2 

0.41991 

RDIF 

1 

NSN 

11 

MA  3 

0 . 27372 

(A) 

MA  6 

0.21114 

MEAN 

8. 1945 

NSN 

12 

AR  3 

0 . 16814 

(A) 

AR  6 

-0. 34273 

MEAN 

26 . 19 

( A) 

=  MUKLUKS 

(B)  -  SOCKS 

(C)  =  COATS 
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the  except  ion  of  a  seasonal  parameter  at  lag  6  for  NSN  5. 
This  means  that  the  seasonality  of  the  data  was  less  import 
ant  in  modeling  the  data  than  other  factors.  Table  VIII 
lists  the  models  in  their  mathematical  notation. 

As  in  the  case  of  the  exponential  smoothing  models  an 
attempt  was  made  to  find  similar  models  between  like  items 
This  failed.  The  demand  patterns  of  like  items  were  not 
similar  enough.  The  models  for  mukluk.  demand  data  include 
differenced  autoregressive  1  and  2,  and  two  models  with 
terms  at  3  and  6  periods,  the  first  in  moving  averages,  the 
second  in  autoregressive  terms.  The  sock  demand  data  model 
are  less  different  in  that  both  include  moving  average 
terms,  although  one  has  only  a  single  moving  average  term 
and  is  undifferenced  and  the  other  has  two  moving  average 
terms  and  required  differencing.  Only  one  model  for  coat 
demand  data  was  generated,  so  no  comparisons  between  coat 
models  could  be  made. 

Because  of  the  increased  complexity  of  the  AP.IMA 
models,  the  models  are  less  likely  (as  was  seen  in  the  case 
of  exponential  smoothing)  to  be  transferable  from  one  item 
to  another  in  spite  of  similar  demands.  The  more  closely 
the  demand  history  is  fitted  by  the  equations,  the  less 
likely  that  equation  will  be  adequate  or  even  useable  on 
another  item,  as  was  seen  with  the  exponential  smoothing 
mode  1 s  . 


Table  VIII.  Mathematical  Notation  for  Box-Jenkins  Models 


NSN  1  (1-01B  -  d-jB^  )  ( 1-B)  Yfc  =  et 

Yt  "  '469Yt-l  -  -167  Yt-2  +  -364  Yt-3  +  et 

MSN  3  Y  -  (  I-6,B)et  +  ii 

Yt  =  et  -  . 289et  +  104.38 

NSN  5  [ll-dQB3  -  d-B6) (1-B) ] Y,  =  ef(l-e,B  -0->B-,) 

J  O  L  L  -L  4.  4. 

Yt  =  Yt_1  -  -624Yt_3  -  . 624Yt_4  +  .478Yt_6 
-  .478Yt_?  +  et  -  •427et_1  +  ,263et_2 

NSN  9  f  ( l-«5 2 B  -  d3B3)]Yt  =  e^l-e^.-,)  +  u 

Yt  =  .  390Yt_1  +  ,261Yt_3  +  et  +  .  246<?t_-. 

+  20.04 


NSN  10  (l-BlYt  -  et(l-eiB  -  ©2B2) 

Y„  -  Yk  t  e.  -  .652e„_1  -  .420e._.-, 


t-1 


t-1 


NSN  11  Yt  =  et ( 1-©3B  -  BgBg)  +  u 

Yt  -  et  -  -274et_3  -  -211et.6  +  8.19 

NSN  12  ((l-d3B3  -  d6B6)]Yt  =  ef  +  u 

Yt  =  ,168Yt_3  -  . 343Yt_g  +  efc  +  u 


Forecasts 

The  forecasts  from  the  various  forecasting  techniques 
are  presented  in  Attachments  9.  10,  11.  and  12  for  December 
data,  January  data.  February  data  and  March  data  respective¬ 
ly.  The  December  data  was  used  to  forecast  demands 
for  January  through  April,  the  January  data  was  used  to 
forecast  demands  for  February  through  April,  and  so  on. 

Mote  that  by  their  very  nature,  the  simple  exponential 
smoothing  and  SBSS  methods  make  a  single  forecast  which  is 
projected  forward  for  all  months.  The  Winters'  and  ARIMA 
models,  however,  project  forward  different  monthly  forecasts 
due  the  their  explicit  consideration  of  seasonality.  At¬ 
tachment  14  contains  graphs  of  the  forecasts  plotted  against 
each  other  and  the  holdout  sample. 

Model  Analysis 

The  different  models,  simple  exponential  smoothing. 
Winters'  seasonal  exponential  smoothing,  ARIMA.  and  SPSS . 
were  evaluated,  as  was  mentioned  in  Chapter  III.  by  how  many 
times  each  model  had  the  lowest  MSE .  In  other  words,  for 
each  NSN.  the  number  of  times  each  model  outperformed  the 
others  for  each  of  the  four  forecast  periods  (December. 
January.  February,  and  March)  was  determined.  The  number  of 
times  each  model  had  the  best  (lowest)  MSE  is  recorded  in 
Table  IX.  The  MSE  values  used  to  determine  the  values  in 
Table  IX  are  listed  in  Table  X  through  Table  XIII.  Recall 
that  only  sev»n  items  had  demand  data  which  could  be  used  to 
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Table  IX.  Number  of  times  each  model  was  "best 


EXPONENTIAL  SMOOTHING 

SIMPLE  WINTERS' 

BOX— JENKINS  SBSS 

NSN  1 

0 

1 

i 

MSN  3 

0 

4 

0 

0 

NSN  5 

0 

0 

4 

0 

NSN  9 

1 

-> 

0 

1 

NSN  10 

0 

4 

0 

0 

NSN  11 

0 

1 

n 

1 

NSN  12 

0 

0 

2 

2 

TOTALS  1  13  9  5 


make  ARIMA  models,  so  NSNs  2.  4.  6.  7.  and  9  are  not  con¬ 
sidered  in  this  analysis. 

Winters'  seasonal  exponential  smoothing  is  clearly  the 
best  forecast  method  using  this  criterion,  having  the  lowest 
MSE  in  13  of  23  forecasts.  APIMA  models  were  the  seconi 


Table 

X. 

Mode  1  MSE ' s 

far  December 

Forecast 

MSE 

FOR  THE  VARIOUS  MODELS 

DEC  FORECAST 

BOX 

EXPONENTIAL 

SMOOTHING 

WINTERS ’ 

SIMPLE 

SBSS 

JENKINS 

NSN 

1 

23 . 38 

48 . 93 

23.49 

2  2 . 36 

NSN 

3 

2880 . 88 

195 .40 

1653 . 97 

3725 . 54 

NSN 

ir 

13.06 

840 .41 

294 . 39 

276.70 

NSN 

9 

270.41 

268 . 89 

431 . 60 

4  30 . 79 

MSN 

10 

4105 . 9? 

480.75 

831 . 73 

1108.13 

NSN 

11 

16 . 23 

32.51 

22.04 

14.13 

NSN 

12 

198 . 88 

1 7  3 . 7  8 

97 . 04 

46.63 

Table  XI.  Model  MSE ' s  for  January  Forecast 


MSE 

BOX 

JENKINS 

FOR  THE  VARIOUS 
JAN  FORECAST 
EXPONENTIAL 

MODELS 

SMOOTHING 

SESS 

WINTERS • 

SIMPLE 

NSN  1 

7.6-? 

29 . 66 

19.88 

18  . 

3-7 

MSN  3 

35  38.0-7 

341.71 

2257 . 15 

4224  . 

9  3 

NSN  5 

14 . 04 

1248 . 65 

306 . 24 

84 

NSN  9 

367 . 33 

244 . 06 

1527.47 

122  . 

36 

NSN  10 

5339 . 82 

1040 .61 

1247 . 37 

1411  . 

J 

NSN  11 

7.58 

49 . 50 

38 . 00 

17. 

36 

NSN  12 

150 . 64 

234.36 

58.46 

58  . 

36 

Table  XII. 

Model  MS 

!E‘s  for  February 

■  Forecast 

MSE 

FOR  THE  VARIOUS 

MODELS 

FEB  FORECAST 

EXPONENTIAL  SMOOTHING 

BOX 

JENKINS 

WINTERS ' 

SIMPLE 

SESS 

NSN  1 

5 . 50 

4.19 

28.80 

33.5-7 

NSN  3 

4171 . 21 

38.14 

2485 . 59 

6045 . 36 

NSN  5 

1  3 . 68 

82.85 

276.81 

294.94 

NSN  9 

377 . 85 

46.03 

111.34 

201.23 

NSN  10 

5127.95 

604 . 16 

1170 . 10 

1897.49 

NSN  11 

5 . 66 

17.20 

22.81 

18 . 87 

NSN  12 

81 .47 

87.87 

92 . 90 

81 . 88 

Table  XIII. 

Model  MSE ' s  for  March 

Forecast 

MSE 

FOR  THE  VARIOUS 

:  MODELS 

MAP.  FORECAST 

EXPONENTIAL 

SMOOTHING 

BOX 

JENKINS 

WINTERS ' 

SIMPLE 

SBSS 

NSN  1 

11.83 

0 . 86 

46.24 

10.54 

NSN  3 

4399 . 67 

60 . 06 

3058 . 09 

3747 . 70 

NSN  5 

0.00 

29.38 

237. 16 

278 . 39 

NSN  9 

501 . 31 

5 . 38 

4 . 97 

192 . 18 

NSN  10 

4919.62 

340 . 40 

1431 . 1 1 

1268 .53 

NSN  11 

8 . 88 

7 . 95 

16.73 

20 . 68 

NSN  12 

0  .  37 

54.76 

259.53 

88 . 05 
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best,  .at  forecasting  with  the  lowest  MSE  in  9  cases.  SBSS 
had  the  "best"  MSE  5  times  and  simple  exponential  smoothing 
had  the  "best"  MSE  1  time. 

I*  should  be  not-1;!  that  in  a  number  of  cases  the  dif¬ 
ferences  between  MSEs  fir  different  models  were  slight.  For 
example,  fo'r  the  December  forecasts  of  MSN  i  SBSS  had  the 
lowest  MSE  at  22. 36.  ARIMA  was  next  best  with  23.39.  simple 
exponential  smoothing  was  third  best  at  23.49.  and  Winters' 
seasonal  exponential  smoothing  had  the  worst  forecast  with 
an  MSE  of  48.93.  The  first  three  models  had  a  spread  of 
only  1.13.  roughly  52$  of  the  lowest  value.  22.36.  This 
shows  that  in  some  cases  several  models  had  very  similar 
performances . 

Several  comments  about  the  forecasts  for  the  SBSS  and 
simple  exponential  smoothing  models  are  in  order .  The  SBSS 
was  the  best  model  only  for  forecasts  from  December  and 
January  when  demands  were  higher  and  closer  t  :<  the  mean  'if 
all  the  demands  and  more  periods  were  being  forecast.  The 
simpler  models  performed  better  when  they  could  forecast 
close  to  the  mean  of  the  entire  series.  Also,  simple  ex¬ 
ponential  smoothing  had  nearly  the  same  forecasts  as  th* 

SBSS  method  in  many,  although  not  all  cases  (NSN  9  is  a  good 
example  of  the  latter  case).  This  is  consistent  with  the 
notion  that  th  the  SBSS  and  simple  exponential  smoothing 
are  simply  modified  moving  averages. 

Considering  the  difficulty  with  which  ARIMA  models  are 
made,  the  methodology  performed  relatively  poorly.  In  only 
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9  of  28  forecasts  was  it  the  best  method.  In  comparison, 
the  simplest  models  (SBSS  and  simple  exponential  smoothing) 
were  "best"  on  6  of  the  28  forecasts. 

Of  the  9  forecasts  in  which  ARIMA  models  excelled,  four 
were  the  forecasts  for  N8N  5.  This  model,  with  AR  ?•  and  •:> 
terms.  MA  1  and  2  terms,  and  a  regular  differencing  term. 

•was  one  of  the  most  complex  developed  in  this  study.  In 
contrast,  the  simplest  ARIMA  model,  an  MA  1  developed  for 
MSN  3.  was  worse  than  the  Winter?1  model  m  all  four  fore¬ 
casts  for  that  NSN.  Winters'  seasonal  exponential  smoothing 
models  also  generated  the  best  forecasts  all  four  times  f :r 
NSN  10.  and  again,  the  ARIMA  model  was  relatively  simple,  an 
MA  1  and  2.  In  contrast .  for  NSN  9.  Winters'  was  the  best 
m  two  forecasts.  SBSS  and  simple  exponential  smoothing  were 
each  best  once,  and  ARIMA  was  never  the  best  model,  m  spit- 
of  having  two  AR  terms  and  one  MA  term.  For  NSNs  11  and  1  '• . 
the  ARIMA  model  was  the  best  model  in  two  of  the  forecasts. 
Winters'  and  SBSS  were  each  best  for  one  forecast  of  NSN  II. 
while  SBSS  was  best  for  two  forecasts  of  NSN  12.  Finally, 
for  NSN  1.  Winters'  was  the  best  model  in  two  forecasts, 
while  ARIMA  was  best  once  and  SBSS  was  best  once.  For  these 
last  3  NSNs.  the  ARIMA  models  had  either  two  MA  or  two  AR 
terms.  So.  it  appears  that  the  more  complex  ARIMA  models 
better  forecast  the  hold-out  data. 

Several  observations  rnay  be  made  concerning  the  rela¬ 
tive  performance  of  the  models.  First,  several  of  th«-  MSFn 
were  quite  close,  such  as  the  December  forecast  for  NSN  11. 
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So.  ARIMA  did  not  perform  as  poorly  as  it  might  appear  at 
first  glance.  Second,  when  ARIMA  was  well-fitted,  as  it 
apparently  was  in  the  case  of  NSN  5.  the  MSE  was  lower  by 
one  order  of  magnitude  as  compared  to  the  other  models,  so 
when  it  was  able  to  capture  the  demand  data  well,  it  ex¬ 
celled.  Third,  it  must  be  remembered  that  the  3ox-Jenkinr 
methodology  works  best  with  extensive  data  f 60  months  or 
better  yet.  120  months).  In  Chapter  III  it  was  noted  that 
the  models  were  built  with  a  bare  minimum  of  data,  which  was 
artificially  made  continuous  from  two  separate  groups  of 
data.  Fourth,  the  exponential  smoothing  models  used  a 
slightly  different  set  of  data  which  was  continuous.  This 
too  may  have  affected  the  results.  Fifth,  the  demand  pat¬ 
terns  were  not  clearly  seasonal,  and  with  possible  man¬ 
agement  influences  changing  demands  from  year  to  year,  any 
model  may  have  difficulty  generating  reasonably  accurate 
forecasts . 

Considering  these  factors,  it  is  not  surprising  that 
the  results  were  mixed.  One  fact  is  clear,  however. 

Winters'  and  ARIMA  models  do  a  better  job  of  forecasting 
demands  that  follow  somewhat  a  seasonal  demand  pattern  than 
do  the  SBSS  and  simple  exponential  smoothing  models  which 
tend  to  best  forecast  slowly  changing  constant  levels  of 
demand . 
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Cone  1 us ions 

This  chapter  discussed  the  analysis  of  the  data  for 
seasonality  and  the  results  of  the  model  building.  The 
models  were  introduced,  along  with  an  analysis  of  whether  or 
not  like  items  had  similar  parameters  in  the  simple  and 
seasonal  exponential  models.  Additionally  the  forecasts 
generated  by  the  models  were  presented,  along  with  an  analy¬ 
sis  of  how  each  technique  performed  compared  to  the  others 

The  seasonality  tests  gave  mixed  results.  Graphically. 
NSNs  1.  2.  3,  5.  9.  10.  11.  and  12  showed  some  consistency 
of  pattern  with  dual  peak  demands  each  year,  while  MSN  A 
showed  a  plateau  shaped  high  demand  period.  NSNs  6.  7,  and 
p  showed  no  discernable  demand  pattern. 

The  exponential  smoothing  models  were  evaluated  fn 
determine  if  one  model  could  be  used  to  forecast  several 
like  items.  One  could  use  a  single  simple  exponential 
smoothing  model  to  generate  forecasts  for  the  three  sizes  of 
mukluks  and  another  model  to  generate  forecasts  for  the 
three  sizes  of  coats.  In  no  cases  could  a  single  Winters' 
seasonal  exponential  smoothing  model  be  used  to  forecast 
demands  for  different  sizes  of  the  same  type  item. 

The  evaluation  of  which  model  best  forecast  the  hold¬ 
out  data  was  also  presented.  Winters'  seasonal  exponential 
smoothing  proved  better  than  the  ARIMA  models.  However, 
based  on  the  small  sample  size,  the  data  difficulties,  and 
the  closeness  of  the  results  this  is  certainly  not  a 
definitive  conclusion. 
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V .  Conclusions  and  Recommendations 


Introduct 1  on 

This  chapter  discusses  hew  the  results  answer  the 
research  questions  and  presents  some  possible  recommenda¬ 
tions  for  both  further  study  and  action.  Briefly,  the  stud 
determined  that  Winters’  seasonal  exponential  smoothing  and 
ARIMA  models  forecast  demands  for  the  sample  items  better 
than  the  SBS5  method  of  forecasting  demands.  Additionally, 
it  was  found  that  the  Winters'  seasonal  exponential  smooth¬ 
ing  models  gave  better  results  than  the  ARIMA  models. 

Responses  to  Research  Questions 

Research  Question  1.  How  does  the  BBSS  currently  ad¬ 
dress  consumable  seasonal  demand  items?  The  BBSS  does  not 
distinguish  between  items  with  seasonal  demand  patterns  arid 
those  without.  Rather,  it  uses  a  modified  simple  exponen¬ 
tial  smoothing  technique  to  calculate  future  demands  for  al 
consumable  items.  Since  simple  exponential  smoothing  tech¬ 
niques  tend  to  lag  changing  trends  in  demand  and  offer 
varying  degrees  of  responsiveness  to  seasonal  demand  pat¬ 
terns  they  are  not  the  best  suited  for  seasonal  demand 
patterns . 

Research  Question  2.  What  are  alternative  methods  of 
addressing  consumable  seasonal  demand  items  in  the  SBBB ? 

Two  primary  alternative  methods  of  addressing  consumable 
seasonal  demand  items  in  the  BBSS  were  studied:  Winters’ 
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seasonal  exponential  smoothing  models  and  the  ARIMA  models. 
Both  forecast  techniques  are  known  for  their  ability  to 
capture  patterns  missed  by  simpler  moving  average  and  simpl 
exponential  smoothing  techniques.  However,  an  argument  for 
the  use  of  either  seasonal  forecast  model  in  the  BBSS  would 
be  tempered  by  added  complexity  of  either  seasonal  model 
with  its  associated  data  storage  requirements. 

Research  Question  3.  Do  items  currently  identified  as 
seasonal  by  personnel  at  the  base-level  actually  display 
seasonal  tendencies?  Two  responses  to  this  question  can  be 
made  depending  on  the  criteria  used  to  answer  it.  A  graph i 
cal  analysis  showed  the  items  identified  as  seasonal  by 
base-level  personnel  have  some,  although  not  strong,  season, 
al  demand  characteristics.  All  twelve  stock  numbers  studio 
had  definitely  lower  summer  demands  compared  to  winter 
demands.  However,  the  timing  of  peak  annual  demands  and 
demand  curve  shapes  differed  from  one  year  to  another, 
Depending  on  the  item,  peak  demand  might  occur  in  January 
one  year,  followed  by  relatively  steady,  relatively  high  de 
mands.  The  following  year,  possibly  because  of  weather, 
differing  management  policies,  or  some  other  reason,  the 
peak  demand  might  b*  one  month  later,  with  few  demands  m 
the  surrounding  months.  This  weakens  the  argument  that 
demands  are  truly  seasonal. 

The  other  technique  determining  the  presence  of  season¬ 
ality  was  to  use  the  autocorrelation  function  and  determine 
if  an  items  monthly  demands  from  twelve  months  apart  were 


correlated.  The  results  of  this  test  were  essentially 
negative.  Although,  as  was  mentioned  in  Chapter  4.  some  of 
the  items  were  close  to  having  significant  correlations 
between  demands  one  year  apart,  the  demand  data  was  not 
clearly  seasonal  according  to  this  test.  A  number  of  rea¬ 
sons  can  be  given  to  account  for  this.  Large  organizational 
buys,  management  control  of  the  items,  the  limited  demand 
history  available  for  study,  fluctuating  weather,  institu¬ 
tional  changes  in  requirements,  low  demand  rates,  and  chang¬ 
ing  budgets  all  could  have  affected  demands  in  a  non -random 
fashion . 

Research  Question  4.  Can  demand  for  the  items  current  - 
ly  identified  as  seasonal  by  personnel  at  the  base-level  be 
better  forecast  using  a  model  that  incorporates  seasonal 
demand  pattern  information,  or  the  current  SBSS  model?  Even 
with  items  that  do  not  show  strong  seasonal  patterns,  using 
forecast  methods  which  can  compensate  for,  or  incorporate 
seasonal  demand  pattern  information,  did  improve  the  ac¬ 
curacy  of  forecasts  as  measured  by  MSE .  Of  28  forecasts  ^f 
7  different  NSNs  the  seasonal  models  (Winters'  and  ARIMA) 
had  the  lower  MSE  22  times.  In  some  cases  there  was  little 
difference  in  the  forecasts;  in  others  the  best  model  per¬ 
formed  much  better.  In  light  of  the  surprisingly  poor 
performance  of  the  ARIMA  models,  and  the  difficulty  of 
analyzing  the  data  and  creating  the  models.  Winters'  was  a 
better  model  choice.  The  Winters'  seasonal  exponential 
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smoothing  models  were  the  best  in  13  of  the  28  forecasts  and 
much  easier  to  develop. 

Research  Question  5.  What  method  could  be  used  by  per¬ 
sonnel  at  the  base-level  to  identify  items  as  seasonal  and 
forecast  demands  for  those  items  accordingly?  Given 
sufficient  demand  data,  a  check  of  the  autocorrelation 
function  is  relatively  easily  accomplished  with  a  microcom¬ 
puter  program  and  can  give  consistent  results.  However, 
given  current  demand  data  storage  requirements  (essentially 
one  year  at  the  base  level)  any  type  of  analysis  of  season¬ 
ality  would  be  impossible  until  several  (4  at  an  absolute 
minimum)  years  of  demand  data  were  collected  for  analysis. 

It  would  be  possible  to  write  into  the  supply  computer  a 
program  that  could  automatically  calculate  the  autocorrela¬ 
tion  function  and  check  at  the  11th.  12th.  and  13th  lags  for 
evidence  of  seasonality,  and  flag  the  item  as  seasonal, 
given  sufficient  demand  data.  To  manually  select  items  for 
review  to  determine  if  seasonality  existed  in  the  demand 
data  would  be  impossible  given  the  number  of  items  in  the 
average  base  supply  account. 

Recommendat ions 

A  number  of  possibilities  exist  for  additional  analysis 
and  study.  In  addition,  some  general  recommendations  are 
made  . 

1.  Recommend  a  follow-on  study  be  accomplished  with 
additional  historical  data.  This  could  answer  many 
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questions  about  the  effect  of  the  small  data  sample  avail¬ 
able  and  the  effect  of  using  non-sequent ia 1  data  in  the 
ARIMA  models.  A  comparison  of  both  models  and  forecasting 
ability  with  additional  data  could  prove  informative. 

2.  Recommend  a  random  sample  of  stock  numbers,  not 
previously  identified  as  seasonal,  be  studied  to  determine 
if  seasonal  techniques  would  be  useful  in  forecasting  de¬ 
mands  and  if  items  not  flagged  as  seasonal  exhibit  seasonal 
demand  patterns.  This  would  be  facilitated  by  the  develop¬ 
ment  of  a  Fortran  computer  program  that  would  take  transac¬ 
tion  history  data  tapes  and  have  as  a  product  useable  month¬ 
ly  demand  data  and  a  monthly  5BSS  forecast  for  comparison 
purposes . 

3.  Recommend  AFLMC/LGS  study  the  feasibility  of 
storing  demand  data  in  a  more  accessible  format  and  for  a 
longer  duration  in  the  SBSS.  Data  is  much  more  likely  to  be 
used  when  it  is  accessible.  As  transaction  histories  are 
accumulated,  samples  of  demand  data  fitting  various  crit¬ 
eria  could  be  extracted  and  analysed  on  mi crocomputer- bas^d 
programs  at  base-level  accounts. 

4.  Recommend  additional  research  be  done  in  the  area 
of  using  Winters'  exponential  smoothing  or  another  seasonal 
model  which  could  better  forecast  seasonal  demands  on  con¬ 
sumable  items  due  to  end-of-year  funding  constraints,  annual 
weather  effects  on  aircraft  parts,  and  annual  exercise 
patterns . 
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5.  Recommend  future  studies  between  the  relationship 
of  more  complex  ARIMA  models  and  SBSS  demand  data  be 
evaluated  to  determine  if  one  to  three  MA  or  AR  term  models 
are  less  capable  of  forecasting  demand  than  Winters'  sea¬ 
sonal  exponential  smoothing. 

.Conclusions 

This  study,  in  spite  of  the  small  sample  size,  and 
possible  data  problems,  demonstrated  that  seasonal  demand 
items  are  better  forecast  using  models  which  account  for 
seasonal  demands.  Although  the  data  requirements  for 
creating  and  using  seasonal  models  are  more  extensive  than 
for  simple  exponential  smoothing  or  moving  average  models, 
with  the  current  reduced  costs  of  storing  data  and  increased 
computer  capabilities,  the  use  of  seasonal  models  is  well 
worth  further  investigation.  The  potential  cost  savings 
resulting  from  better  forecasts  of  consumable  items  certain¬ 
ly  warrant  further  investigations  into  this  area. 
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Attachment  1:  Items  Selected  for  Study 


STOCK  NUMBER 

NOMENCLATURE 

LIKE  ITEM 
GROUP 

NSN1 

8430002690100 

BOOT.  MUKLUK  LARGE 

(A) 

NSN2 

8415003761661 

PARKA. ECW  GRN  X-SMA 

NSN3 

8440001536717 

SOCKS. MENS  C/W  NA10 

(B) 

NSN4 

8415007822935 

COAT  MANS  SML-SHOR 

(C) 

NSN5 

8415007822938 

COAT  MAN  MED  SHT 

(C) 

NSN6 

8440001536721 

SOCKS', MENS  C/W  NA14 

(B) 

NSN7 

8415007822940 

COAT  MAN  MED  LG  33 

(C) 

NSN8 

8415002702014 

UNDERSHIRT  EX  WE A  L 

NSN9 

8415002687872 

GLOVE  SHELLS.  CW  SZ2 

NSN10 

8440002614897 

SOCKS . MEN , CW . NAUTR9 

(B) 

NSN1 1 

8430002690098 

BOOT  MUKLUK  MANS  SM 

(A) 

MSN  12 

8430002690099 

BOOT  MUKLUK  MED 

( A ) 
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Attachment  2:  Sample  Data  Extraction  Fortran  Program 

and  Data  Tape  Schema 


DEFINE  VARIABLES 

CHARACTER* 2  TTPC 

CHARACTER* 3  DIG 

CHARACTER* 4  DOLD . DOLT 

CHARACTER  *  6  TRAN3D . ACTQTY 

CHARACTER  *  1 3  3TKNUM . SNREQ . STOCKN 

CHARACTER* 14  DOCNUM 

CHARACTER *20  FILENM 

CHARACTER* 121  DUMMY 

INTEGER *4  I COUNT . OCOUNT 


SET  RECORD  COUNTERS 

I COUNT- 0 
OCOUNT-O 

GET  INPUT  FILE  NAME  AND  OPEN  FILE 
5  WRITE (*.900) 

900  FORMAT ('  ENTER  INPUT  FILE  NAME:  '.$) 

READ (*,905)  FILENM 
905  FORMAT (A20) 

OPEN ( 1 , FILE=FILENM . STATUS = ' OLD ' ) 

GET  OUTPUT  FILE  NAME  AND  OPEN  FILE 
10  WRITE (*.910) 

910  FORMAT ( 1  ENTER  OUTPUT  FILE  NAME:  '.$) 

READ (*.915)  FILENM 
915  FORMAT (A20) 

GET  STOCK  NUMBER  FOR  SEARCH 

WRITE ( * .917) 

917  FORMAT ( 1  ENTER  STOCK  NUMBER:  ‘.$) 

READ (*.918)  STOCKN 

918  FORMAT ( A13 ) 

READ  A  RECORD 

OPEN ( 2 . FILE=FILENM . STATUS* 'MEW  ) 

2  0  READ (1.920. END- 9 9  9 ) 3TKNUM , DUMMY . D I C . DUMMY . 
&  DOCNUM . DOLD . DUMMY . TRANSD . DUMMY . 

&  ACTQTY . DUMMY . DOLT . DUMMY . TTPC . 

&  DUMMY . SNREQ 

I COUNT = I COUNT + 1 


1 14- 


920  FORMAT ( A13. A13. A3.A13. A14. A4  . A6.A6.A9. 

&  A6.A8.A4.A8.A2.A16.A13) 

C 

C  STOCK  NUMBER  SELECTION  LIST 
C 

I F ( STKNUM . EQ . STOCKN )  THEN 

WR ITE ( 2 . 940 ) STKNUM . D I C . DOCNUM . DOLD . 

*  TRANS D . ACTQTY . DOLT . TTPC . SNREQ 

940  FORMAT (1X.A13.2X.A3.2X.A14.2X.A4.2X. 

&  A6.2X.A6.2X.A4.2X.A2. 1X.A13) 

OCOUNT-OCOUNT+1 
END  IF 
r 

GOTO  20 
C 

999  WRITE (2.950) I COUNT . OCOUNT 

950  FORMAT ('  INPUT  RECORDS  READ  =  '.15. 

&  '  OUTPUT  RECORDS  WRITTEN  =  .15) 

CLOSE ( 1 ) 

CLOSE ( 2 ) 

STOP 

END 
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Data  Tape  Schema  (27) 


POS  . 

SECTOR 

LA 

DESCRIPTION 

NO. 

POS. 

15 

1-15 

SN 

Stock  Number 

2 

16-17 

WS 

System  Designate! 

i 

18 

AC 

Srar.  Type 

3 

19-21 

ER 

ERR  CD 

1 

22 

ME 

Stockage  Priority  Code 

2 

23-24 

PY 

Issue  Priority 

1 

25 

TX 

TEX 

1 

26 

PM 

Demand  Code 

3 

27-29 

ID 

DIC-TRIC 

2 

30-31 

UI 

Unit-of-Issue 

2 

32-33 

FC 

Fund  Code 

6 

34-39 

SA 

Sup-Requ is i t loner 

3 

40-42 

RI 

Routing  Identifier 

14 

43-56 

DN 

Document  Number 

4 

57-60 

LD 

DOLD 

6 

6  i-66 

EB 

Ending  Balance 

6 

67-72 

DT 

Transaction  Date 

6 

73-78 

TN 

Transaction  Ser  Nbr 

3 

79-81 

FI 

FIA  Trans 

6 

82-87 

QY 

Action  Qty 

8 

88-95 

XP 

Extended  Cost 

4 

96-99 

LT 

DOLT 

3 

100-102 

ST 

Status  or  Advice  Code 

1 

103 

Blank 

3 

104-106 

OP 

Output  Terminal 

1 

107 

FO 

MAT.  CAT.  SOS  Code 

2 

108-109 

PE 

TTPC 

1 

110 

IN 

Print  Punch.  Old  TR 

111 

BJ 

Budget  Code 

14 

112-125 

MF 

Mark-For 

15 

126-140 

SR 

SN  Requested 

19 

141-159 

NM 

Noun 

5 

160-164 

MC 

MFC ’ s  Code 

1 

165 

RW 

Reason-Why-Code 

7 

166-172 

Blank 

i 

173 

RC 

Reporting  Code 

1 

174 

IX 

LEX 

19 

175-193 

Blank 
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DATE 

JAN  82 

NS'Nl 

MSN  2 

NSN3 

NSN4 

NSN5 

NSN6 

FEB 

0 

0 

0 

0 

0 

0 

MAR 

4* 

0 

48 

20 

23 

0 

APR 

3 

0 

2  3 

o  9 

28 

0 

MAY 

0 

0 

0 

0 

0 

0 

JUN 

0 

0 

20 

4 

4 

0 

JUL 

0 

0 

9 

1 

7 

0 

AUG 

8 

0 

90 

5 

7 

0 

SEP 

0 

0 

197 

18 

9  9 

14 

OCT 

13 

4 

1040 

83 

66 

100 

NOV 

104 

0 

206 

69 

0 

13 

DEC 

0 

0 

0 

0 

0 

0 

JAN  83 

101 

4 

267 

23 

68 

0 

FEB 

MAR 

120 

3 

112 

42 

8 

0 

APR 

0 

0 

186 

10 

12 

0 

MAY 

0 

0 

8 

4 

5 

5 

JUN 

0 

0 

15 

9 

7 

0 

JUL 

AUG 

0 

0 

2 

2 

4 

0 

SEP 

0 

0 

160 

14 

10 

0 

OCT 

0 

0 

0 

13 

26 

] 

NOV 

0 

0 

47 

5 

12 

~) 

J 

DEC 

JAN  84 

0 

3 

20 

26 

17 

14 

FEB 

MAR 

APR 

MAY 

JUN 

JUL 

AUG 

SEP 

OCT 

NOV 

DEC 

JAN  8 “5 

FEB 

MAR 

APR 

MAY 

JUN 

JUL 

AUG 

SEP 

OCT 

NOV 

DEC 


117 


DATE 

JAN  86 
FEB 

MAR 

APR 

MAY 

JUN 

JUL 

MSN  1 

NSN2 

AUG 

3 

1 

SEP 

55 

1 

OCT 

NOV 

DEC 

JAN  87 

6 

6 

FEB 

MAR 

37 

8 

APR 

17 

6 

MAY 

1 

0 

JUN 

0 

0 

JUL 

0 

0 

AUG 

5 

0 

SEP 

2 

0 

OCT 

32 

9 

NOV 

11 

3 

DEC 

JAN  88 

4 

0 

FEB 

47 

2 

MAR 

6 

3 

APR 

0 

1 

MAY 

0 

0 

JUN 

0 

0 

JUL 

0 

0 

AUG 

1 

0 

SEP 

6 

3 

OCT 

27 

Q 

NOV 

4 

0 

DEC 

8 

1 

JAN  89 

7 

0 

FEB 

10 

1 

MAR 

9 

0 

APR 

5 

1 

NSN3  NSN4  NSN5  NSN6 


o 

9 

1 

0 

141 

30 

31 

0 

146 

13 

31 

2 

211 

44 

58 

21 

42 

22 

21 

4 

80 

3 

9 

0 

4 

0 

i 

0 

3 

2 

3 

n 

98 

8 

5 

0 

9 

3 

3 

0 

325 

19 

16 

93 

533 

4 

107 

20 

0 

16 

7 

0 

32 

18 

139 

y 

38 

12 

4 

9 

13 

12 

0 

0 

22 

6 

4 

0 

4 

0 

1 

o 

0 

1 

1 

3 

40 

5 

3 

14 

236 

13 

10 

57 

154 

13 

6 

o 

74 

1 

13 

■*> 

73 

8 

6 

29 

77 

5 

5 

0 

63 

6 

1 

7 

49 

1 

0 

Q 

25 

1 

1 

0 
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NSN7 

NSN8 

NSN9 

NSN10 

NSN11 

NSN12 

DATE 

JAN  82 

FEB 

0 

0 

0 

0 

0 

0 

MAR 

24 

111 

30 

45 

10 

13 

APR 

10 

76 

8 

7 

9 

MAY 

0 

n 

0 

0 

0 

0 

JUN 

3 

5 

~7 

/ 

10 

6 

6  ' 

JUL 

a 

i 

8 

31 

0 

0 

0 

AUG 

5 

220 

34 

81 

5 

10 

SEP 

28 

216 

24 

207 

1 

0 

OCT 

71 

3010 

29 

310 

n 

/ 

23 

NOV 

20 

355 

53 

159 

54 

105 

DEC 

0 

0 

0 

0 

0 

0 

JAN 

10 

515 

59 

70 

28 

20 

FEB  83 

27 

243 

105 

311 

1 

1  36 

MAR 

6 

106 

7 

41 

1 

0 

MAY 

1 

12 

59 

0 

0 

o 

JUN 

9 

11 

9 

1 

0 

0 

JUL 

0 

5 

3 

18 

0 

0 

AUG 

11 

28 

19 

158 

0 

o 

OCT 

15 

1008 

30 

127 

14 

168 

NOV 

2 

6 

0 

12 

0 

0 

DEC 

14 

249 

18 

138 

1 

51 

JAN  84 

FEB 

MAR 

APR 

MAY 

JUN 

JUL 

AUG 

SEP 

OCT 

NOV 

DEC 

JAN  85 

FEB 

MAR 

APR 

MAY 

JUN 

JUL 

AUG 

SEP 

OCT 

NOV 

DEC 

JAN  86 
FEB 
MAR 
APR 
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NSN7  NSN8  NSN9  NSN10  NSN11  NSN12 


DATE 

MAY 

JUN 

JUL 

AUG 

0 

12 

0 

9 

2 

6 

SEP 

30 

342 

92 

107 

18 

6  6 

OCT 

14 

265 

15 

80 

12 

5  ij 

NOV 

DEC 

JAN  87 
FE3 

18 

357 

12 

129 

37 

79 

MAR 

APR 

3 

49 

14 

61 

O 

La 

9 

MAY 

2 

9 

18 

6 

2 

4 

JUN 

1 

14 

2 

9 

1 

0 

JUL 

0 

91 

6 

6 

1 

2 

AUG 

5 

216 

3 

80 

x 

11 

SEP 

9 

23 

10 

'“7 

5 

3 

OCT 

li 

158 

9 

69 

30 

44 

NOV 

li 

590 

1 

206 

4 

o 

DEC 

i 

86 

23 

78 

0 

n 

JAN  88 
FEB 

95 

27 

2 

18 

0 

63 

MAR 

3 

153 

4 

62 

5 

24 

APR 

9 

La  La 

0 

2 

0 

3 

MAY 

4 

55 

4 

4 

28 

4 

JUN 

1 

4 

0 

1 

JUL 

0 

n 

La 

1 

5 

n 

1 

AUG 

2 

42 

2 

29 

ij 

0 

SEP 

15 

385 

56 

100 

27 

4  2 

OCT 

5 

412 

72 

277 

30 

65 

NOV 

4 

2 

4 

2 

14 

25 

DEC 

3 

502 

13 

17 

6 

11 

JAN  89 

2 

873 

52 

81 

6 

23 

FEB 

1 

161 

12 

33 

4 

18 

MAR 

0 

54 

■y 

28 

4 

31 

APR 

1 

47  ' 

3 

17 

5 

13 
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Attachment  4:  Exponential  Smoothing  Data 
( Under  1 i ne  marks  values  created  using  Delphi  technique) 


DATE 

NSN1 

NSN2 

NSN3 

NSN4 

NSN5 

NSN6 

AUG  86 

3 

1 

9 

9 

1 

0 

SEP 

55 

1 

141 

30 

31 

0 

OCT 

6 

f. 

146 

13 

31 

•"> 

NOV 

18 

3 

124 

') 

0 

DEC 

23 

3 

140 

33 

19 

12 

JAN  87 

39 

3 

153 

36 

28 

13 

FEB 

37 

8 

211 

44 

58 

-i  i 

_  j- 

MAR 

18 

6 

92 

22 

26 

10 

APR 

17 

6 

42 

22 

21 

4 

MAY 

1 

0 

80 

3 

7 

0 

JUN 

0 

0 

4 

0 

1 

0 

JUL 

0 

0 

3 

2 

3 

■1 

AUG 

5 

0 

98 

8 

5 

0 

SEP 

2 

0 

Q 

3 

3 

o 

OCT 

32 

2 

325 

19 

16 

93 

NOV 

11 

3 

533 

4 

107 

20 

DEC 

4 

0 

0 

16 

7 

0 

JAN  88 

24 

64 

16 

39 

x 

FEB 

47 
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Attachment  5:  Box-Jenkins  Data 
(Under  1 i ne  marks  values  created  using  Delphi  technique) 
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Attachment  6:  Box-Jenkins  Forecast 
Demand  Data  Plots 
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Figure  26.  Box-Jenkins  Forecast  Demand  Data  for  NSN 
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gure  28.  Box-denkins  Forecast  Demand  Data  for  NSN 
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Figure  29.  Box-Jenkins  Forecast  Demand  Data  for  NSN 
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Box-Jenklns  Forecast  Demand  Data  for  NSN  10 
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Figure  31.  Box-.Jenkins  Forecast  Demand  Data  for  NSN  11 
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gure  32.  Box-Jenkins  Forecast  Demand  Data  for  N5N  12 


Attachment  7:  Autocorrelation  Function  Plots 
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33.  Autocorrelation  Function  for  NSN 
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gure  34.  Autocorrelation  Function  for  NSN 
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Autocorrelation  Function  for  NSN 
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36.  Autocorrelation  Function  for  NSN 
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Figure  37.  Autocorrelation  Function  for  NSN  6 
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re  38.  Autocorrelat ion  Function  for  NSN 
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Figure  39.  Autocorrelation  Function  for  NSN 
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Figure  40.  Autocorrelation  Function  for  NSN 


uoueidjjoxjyiv 
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Figure  41.  Autocorrelation  Function  for  NSN  10 
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on  Function  for  NSN  11 
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Attachment  8:  ARIMA  Models  Computer  Output 


THE  ESTI KITED  RESIDUES  --  KOOEL  1 
CURUUT1VE  PED10D06RU  1  PIOUIIIITT  UR1TS 

no  .100  150  200  250  300  350  .400  450 


Figure  44.  Cumulative  Penodogram  for  NSN  1 
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Figure  45.  Estimated  Residuals  Histogram  for-  NSN  1 
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THE  esrifWTED  RESIDUALS  --  KflOEL  1 
GRAPH  OF  OBSERVED  SERIES  ACE 
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Figure  46.  Residual  ACF  for  NSN  1 
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Figure  47.  Log  Spectrum  for  NSN  1 
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Figure  48.  C'imulative  Penodogram  for  NSN  3 
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Figure  49.  Estimated  Residuals  Histogram  for  NSN  3 


149 


THE  EST1HATE0  RESIDUALS  --  MODEL  1 
6RAPH  OF  OBSERVED  SERIES  ACF 

5.000  10.000  15.000 


.900 


.750 


.600 


.450 


.300 


.150 


20.000  25.000  30.000  35.000  40.000  45.000 


1.050  ♦ 

£  0 

I 
* 

X 
* 

X 
I 
X 
I 
♦ 

I 
I 
I 
I 
♦ 

I 
X 
I 
I 
♦ 

I 
X 
I 
X 
♦ 

I 
X 
X 
I 

t 

I 
X 
I 
X 

.000  ♦  0  0  0  0  ♦ 


I 


1 1 


I  X 

I I 
I 

X 


10  0  0 
1 1 1 
I  I  X 

I  X  I 

I I 


X 

X  I 
♦  0 


I 

X 

X  0  0  0 
X  I 


I 


1 1 
I  X 

X  XIX 

I  III 

0  *  0  0  0  0  * 
X  X 
X 


-.150  ♦ 


-.300  ♦ 


100*0000 
X  X  X  X  I  I  I 

X  I  X 

I  I  I 


I 

I 

I 

I  I 
♦  0 


I10*0000*0000* 


Figure  50.  Residual  ACF  for  NSN  3 
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Figure  51.  Log  Spectrum  for  NSN  3 
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Figure  52.  Cumulative  Periodogram  for  NSN  5 
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Figure  53.  Estimated  Residuals  Histogram  for  NSN  5 


153 


THE  ESTIMATES  RESIDUALS  --  MODEL  I 
6RAPH  Of  OBSERVED  SERIES  ACE 

5.000  10.000  15.000  20.000  25.000  30.000  35.000  <0.000  <5.000 


1.050  ♦♦♦♦♦♦♦♦♦♦♦ 
E  0 

I 

I 

I 

.900  ♦♦♦♦♦♦♦♦•♦»♦ 
I 
I 
I 
I 

.750  ♦  ♦♦♦♦♦♦♦♦♦♦ 

I 
I 
I 
* 

.600  ♦♦♦♦♦♦♦♦♦♦♦ 
I 
I 
I 
I 

.<50  ♦  ♦♦♦♦♦♦♦♦♦♦ 

I 
I 
X 
X 

.300  ♦  ♦♦♦♦♦♦♦♦♦♦ 

l 
( 

I 

X 

.150  ♦  ♦  ♦  ♦  I  ♦  ♦  ♦  f  ♦  ♦  ♦ 

I  I  1 

I  I  I  I 

II  I  II  I 

I  l  I  I  I  I  I  III 

.OOO*OOOO*IOOO*OOIO*OOOC*OOOOMOOO*OOOOMO0O»OOOO*OOOO* 
III  II  II  III 

I  I  II  III 

II  II  I 

II  II  I 

-.150  ♦  *  ♦  ♦  ♦  I  ♦  ♦  ♦  ♦  ♦  ♦ 

I  I 

I 


-.300  ♦♦♦♦♦♦♦♦♦♦♦ 


Figure  54.  Residual  ACF  for  NSN  5 
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Figure  55.  Log  Spectrum  for  NSN  5 
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Figure  56.  Cumulative  Penodogram  for  NSN  9 
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Figure  57.  Estimated  Residuals  Histogram  for  NSN  9 
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Figure  58.  Residual  ACF  for  NSN  9 
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Figure  59.  Log  Spectrum  for  NSN  9 
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Figure  60.  Cumulative  Periodogram  for  NSN  10 
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Figure  61.  Estimated  Residuals  Histogram  for  NSN  10 
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Figure  63.  Log  Spectrum  for  NSN  10 
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Figure  64.  Cumulative  Periodogram  for  NSN  11 
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Figure  65.  Estimated  Residuals  Histogram  for  NSN  11 
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Figure  69.  Estimated  Residuals  Histogram  for  NSN  12 
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Attachment  10:  December  Forecasts 
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Attachment  10: 
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Forecasts 
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Attachment  11:  February  Forecasts 
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February  Forecasts 
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Attachment  12:  March  Forecasts 
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Attachment  13:  Forecast  Graphs 
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gure  72.  December  Forecasts  for  NSN 
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Figure  73.  December  Forecasts  for  NSN  3 
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74.  December  Forecasts  for  NSN 
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75.  December  Forecasts  for  NSN 
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gure  76.  December  Forecasts  for  NSN  10 
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igure  77.  December  Forecasts  for  NSN  11 
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78.  December  Forecasts  for  NSN  12 
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gure  79.  January  Forecasts  for  NSN 


Model  Forecasts  vs.  Actuals 


■*se:»jOd  puButao 


gure  80.  January  Forecasts  for  NSN 
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gure  81 .  January  Forecasts  for  NSN 
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gure  fi2.  January  Forecasts  for  NSN 


Figure  83.  January  Forecasts  for  NSN  10 
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85.  January  Forecasts  for  NSN  12 
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Figure  86.  February  Forecasts  for  NSN 
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gure  88.  February  Forecasts  for  NSN  5 
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Figure  89.  February  Forecasts  for  NSN  9 
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Figure  90.  February  Forecasts  for  NSN  10 
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gure  91.  February  Forecasts  for  NSN  11 
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gure  92.  February  Forecasts  for  NSN  12 
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